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A FREQUENCY DOMAIN APPROACH TO THE DESIGN 
AND ANALYSIS OF LINEAR MULTIVARIABLE SYSTEMS 

By William A. Wolovich* 

Office of Control Theory and Application 
Electronics Research Center 

SUMMARY 

A new approach to the design and analysis of linear multivari- 
able systems is given. Several applications of one rather simple 
basic result, called the structure theorem, are presented. This 
theorem allows us to combine time and frequency domain information 
in a compact, concise expression for the transfer matrix of a 
linear multivariable system. This expression and the definitions 
required to formulate it are used to establish new results, which, 
for the most part, can be classified as design algorithms that can 
be implemented by using pencil and paper methods alone, if the 
system order is relatively low, or by employing rather simple 
computer programs. 

In particular, a new technique for obtaining realizations, 
given the transfer matrix of a system, is introduced. A direct, 
constructive technique for arbitrarily assigning all eigenvalues 
(closed loop poles) of a controllable multivariable system by us- 
ing linear state variable feedback is also given. A simplified 
expression for the characteristic equation of a system compensated 
by linear output feedback is discussed from the viewpoint of pole 
assignment. Questions pertaining to linear optimal control are 
considered, and improvements are made over prior work in the areas 
of spectral factorization and the root-square locus. Specifically, 
a solution to the output regulator problem is formulated in the 
frequency domain, and spectral factorization is employed to obtain 
the optimal control law solution. A simplified expression is given 
for formulating and plotting the root-square locus. Various ques- 
tions pertaining to the design of noninteracting control systems 
are also considered, and a procedure is given for determining which 
poles cannot be altered while decoupling a system using linear 
feedback. A design algorithm for achieving n aximum pole assign- 
ment under linear decoupling feedback is presented. All systems 
which can be decoupled by using input dynamics are characterized, 
and an algorithm is outlined for achieving an asymptotically 
stable decoupled design with arbitrary pole placement. Two design 
algorithms for achieving a noninteractive design via linear state 
variable feedback are given. The design of a helicopter stability 
augmentation system based on desired handling qualities is pre- 
sented. Many of the results are applied, and pencil and paper 
methods are employed to achieve the final control system used. 


*The research reported has been carried out under the supervision 
of Prof. Peter L. Falb and has also been submitted to Brown Univer- 
sity as a dissertation in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy in the Electrical Sciences. 
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1. INTRODUCTION 


We have witnessed a radical change in the methods used for 
the design and analysis of control systems during the past decade 
Control systems engineers, who once relied almost entirely on 
classical frequency domain techniques such as the root locus. 

Bode plots, and Nyquist diagrams, are gradually adopting the 
modern state space approach to system design. This is primarily 
due to the rapid theoretical advances made over the past few 
years, especially in the area of optimal control, and to the evo- 
lution of the high speed, easily accessible digital computer. 
Designs which once took weeks of effort involving excessive multi 
loop analysis and elaborate trial and error simulation can now be 
obtained in a matter of minutes or even seconds. Today's control 
literature is filled with applications of modern control theory, 
via computational algorithms, to the design and analysis of dy- 
namical systems and the end is not in sight. However, we cannot 
be overly optimistic in evaluating the impact of modern control 
theory on practical system design. Many systems are still too 
complex for even the latest theory or computational algorithms to 
handle, and trial and error simulation and evaluation appear to 
be the only alternative in achieving a satisfactory design. 
Furthermore, computational algorithms employing modern control 
principles do not always lead to a satisfactory design or a de- 
sign which can easily be altered. Frequency domain techniques 
still retain certain advantages over state space methods, espe- 
cially for single input, single output systems. In the case of 
multivariable systems, however, time domain methods have thus far 
been dominant. 

The primary purpose of this report is to present new insight 
into the design and analysis of linear multivariable systems. 

This is accomplished by presenting a number of applications of 
one rather simple basic result, called the structure theorem. 

This theorem, which can easily be stated and established, allows 
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us to combine time domain and frequency domain information in one 
concise expression for the transfer matrix of a linear multi- 
variable system. We use this compact transfer matrix expression 
to establish several new results, which represent the main con- 
tribution of this report. These results, for the most part, can 
be classified as design algorithms. An important feature of many 
of these algorithms is that they can be implemented by using 
pencil and paper methods alone, if the system order is relatively 
low. For higher order systems, appropriate computer algorithms 
can be employed to simplify the design process. 

We introduce the structure theorem in Section 2 . Sections 
3, 4, and 5 are, for the most part, mutually independent, and 
represent a variety of research areas where the structure theorem 
has been successfully applied. Section 6, which deals with a 
practical design problem, uses the results obtained in most of 
the prior sections. Each section contains an introduction and 
several examples to assist the reader in understanding the various 
steps taken. A list of symbols and special notations are included 
as Appendixes A and B respectively. 

2. A STRUCTURE THEOREM FOR LINEAR MULTIVARIABLE SYSTEMS 

The primary purpose of this section is to state and prove a 
simple basic result which underlies the remainder of this report. 
This result, which will be called the structure theorem, enables 
one to express the transfer matrix of any linear, time-invariant 
dynamical system in a concise and compact manner. This expression 
for the transfer matrix is based on a canonical form for the 
dynamical equations which characterize the system. The structure 
theorem then provides a means for directly expressing the transfer 
matrix of the system in terms of the canonical equations of motion. 
The effect of linear state variable feedback on the transfer matrix 
of the system is also clarified. All of these points are discussed 
in reference 1 and will be covered in depth in Sections 2.1, 2.2, 
and 2.3. 
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In Section 2.1, we define the class of systems considered 
and discuss linear state variable feedback and its implications. 
Section 2 . 2 serves to introduce the structure theorem for the case 
of systems which are completely controllable. The general struc- 
ture theorem is presented in Section 2.3. 

2.1 Transfer Matrices and State Variable Feedback 

The class of systems which will be considered in this report 
can be defined in the following manner: 

x = Ax + Bu , y = Cx (1) 

where x is an n vector, called the state; u is an m vector, called 
the input; y is a p vector, called the output; and A, B, and C are 
constant matrices of the appropriate dimension. We assume that 
the m inputs and p outputs are linearly independent, i.e., B and 
C are of full rank.* Furthermore, assuming zero initial condi- 
tions on the state x, we can determine and express the open loop 
transfer matrix, T q (s), relating the Laplace transform of the 
output y(s) , to the Laplace transform of the input, u(s) , of 
Eq. (1) , in terms of the triple {A, B, C} as 

T q (s) = C (si - A) -1 B (2) 

We will be interested in the effect of linear state variable 
feedback on the transfer matrix of the system. In particular, 
suppose that 


u = Fx + Gw ( 3 ) 

where w is an m-dimensional external input, and F and G are con- 
stant matrices of appropriate dimension. Substituting Eq. (3) 


*This condition should hold in most cases. If it does not, we 
merely reduce the dimension of B and/or C until it does. 
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for u into Eq. (1), and solving for y( s ) terms of w(s), we 
obtain the closed loop transfer matrix, 

T FfQ (s) = C (si - A - BF)“ 1 BG (4) 

Paramount to the development of the structure theorem is the 
observation that the transfer matrix of a linear system repre- 
sents an input/output relationship and, consequently, should be 
independent of the choice of "state" for the system. In particu- 
lar, if we consider the effect of altering the state, x, of the 
system (Eq. (1)) via an (nxn) nonsingular similarity transforma- 
tion Q, by setting 


z = Qx (5) 

we obtain an "equivalent" system 

z = QAQ -1 z + QBu , y = CQ _1 z (6) 

or 

Cz (7) 

We now establish an irn- 

P 1 : Tr ansfer matrices of equivalent systems are identical . 
Proof : Simply noting that 

C (si - A) _1 B = CQ _1 Q (si - A) _1 Q“ 1 QB 

= CQ -1 (si - QAQ" 1 )~ 1 QB = C (si - A) _1 B 

(8) 

actually go one step further and establish the 
transfer matrices of equivalent systems after the 
of linear state variable feedback. In particular. 


We can 
equality of 
application 


z = Az + Bu , y . = 

where A = QAQ -1 , B = QB, and C = CQ -1 . 
portant well known proposition (Pi) : 
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if u = Fx + Gw = FQ _1 z + Gw = F z + Gw, where F = FQ -1 , then 

-F,G ( s) = £(sl " A - BF) _1 BG = C(sl - A - BF)BG = T* Q (s) (9) 

We note that Eq. (9) is a relatively simple expression for the 
closed loop transfer matrix of the system in terms of its time 
domain representation. However, the reduction of Eq. (9) to a 
(pxm) matrix of transfer functions requires the inversion of the 

A A A 

(nxn) matrix (si - A - BF) or (si - A - BF) . Furthermore, the 
effect of linear state variable feedback on the pm terms compris- 
ing the closed loop transfer matrix is not at all apparent. The 
structure theorem, which will be introduced in the next section, 
significantly corrects these deficiencies. 

In developing the structure theorem for multivariable systems, 
we will first consider single input (m = 1) systems, sometimes 
referred to as scalar systems. The discussion dealing with this 
class of systems is not intended to be rigorous, since all of the 
statements which will be made here either have been established 
elsewhere or will be established for the more general (multi- 
variable) class of systems in Section 2.2. 

For scalar systems, it is well known (ref. 2) that (1) linear 
state variable feedback (u =. fx + w) affects only the closed loop 
poles of the transfer function, and (2) any closed loop configura- 
tion can be achieved via linear state variable feedback if the 
pair {A,b} is completely controllable. In particular, if the pair 
{A,b} is completely controllable, one can find a similarity trans- 
formation Q, which transforms the pair {A,b} to a controllable 
companion form (ref. 3), after which, the effect of state feed- 
back is immediately apparent. This point can easily be illustrated 
by example. Consider the following completely controllable pair 
{A,b} representing a scalar system whose equations of motion are 
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given by Eq. (1) , with 


10 2 

2-27 b = 

113 

The output y (and hence c) can be chosen arbitrarily and plays no 
important role at this point. Suppose we are required to study 
the effects of linear state variable feedback on the closed loop 
transfer function of this system. We note that the closed loop 
poles (zeros of the denominator of the transfer function) of the 
system are equal to the eigenvalues of (A + bf) (ref. 2) . Further- 
more, the numerator of the transfer function is unaffected by f, 
the feedback vector (ref. 2) . The effect on the closed loop poles 
of the system (the eigenvalues of (A + bf) ) , due to variation of 
f, is not at all apparent in terms of the given pair {A,b}. How- 
ever, suppose we transform the pair {A,b} to controllable compan- 
ion form using the similarity transformation Q, where 

1-11 
0 1-2 
0 0 1 

Techniques for finding a suitable Q are given in references 4, 5, 

6, and 7, for completely controllable scalar systems. In terms 
of the new state z = Qx, we obtain the differential equations 
z = QAQ 1 z + Qbu = Az + bu , where 

0 10 0 

A = 0 0 1 b = 0 

10 2 1 

If one now deals with the new state z and the corresponding 

. - A ^ - 

pair iA,b>, it is clear that linear state variable feedback. 
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yields the closed loop matrix 


u = fz + w, where f — 


A + bf 


0 10 

0 0 1 

f.^+1, ^ 2 ' ^3"*"^ 


Clearly, one can select the feedback vector f to achieve any de- 
sired third (final) row of the companion matrix (A + bf) . Since 
the characteristic polynominal of a companion matrix is expressed 
directly in terms of the elements comprising the final row (ref. 8) 
that is, 

A AA . ^ A 0 ^ A 

|sl - A - bf | = s - (f + 2 ) s x - f 2 s - f x - 1 

/\ /\ /N 

it follows that all three eigenvalues of (A + bf) can be arbi- 

A A A “1 

trarily chosen. Furthermore, since (A + bf) = Q (A + bf)Q , 

/V. 

where f = fQ, the eigenvalues of (A + bf) are the same as those 

A A 

of (A + bf) (ref. 8). Therefore, we conclude that the effect of 
linear state variable feedback on the eigenvalues of (A + bf) , 
actually on the coefficients of the characteristic equation of 
the given scalar system, can easily be understood by considering 
the equivalent companion form system. A similar result holds for 
multivariable systems and, as we will show, forms the basis of 
the structure theorem. 


2 . 2 A Structure Theorem for Controllable Systems 

In the case of multivariable systems, we are interested in 
finding a similarity transformation Q, which transforms the given 
system to a canonical form similar to the companion form for 
scalar systems. As we will show, an appropriate Q always exists 
and can be found if the pair {A,B} is completely controllable by 

__ -j 

first considering the (nxnm) matrix K = [B,AB,...,A n X B] . 

Clearly, K has rank n if the pair {A,B} is completely controllable, 
and it is possible to define a "lexicographic" basis for con- 



sisting of the first n linearly independent columns of K possibly 
reordered (ref. 7) . We let L be the matrix whose columns are the 
elements of the "lexicographic" basis so that 


L 


b , Ab.^ 


a l ^ 

- 1 ' -2 



(10) 


where b^,...,b m are the columns of B. Setting 


k 

d Q = 0 , d^ ~ k = 1 , 2 , . . . ,m (11) 

i=l 

and letting be the d^-th row of L we can see that the 
matrix Q given by 


Q 


&1 

£; 1 A 




& A 

-m- 


m 


-1 


-1 


( 12 ) 


represents a similarity transformation for which Eq. (7) is in 


"companion" form (refs. 7 and 9). 
the transformed state Eq. (7), A 
form 


More precisely, in terms of 

A 

(a. .) is a block matrix of the 
ID 


A = 


tn 

-21 


A i 
-ml 


-2m 


A 

-mm 


(13) 
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with a (o\xcr^) companion matrix given by 



and A. . a 
-iD 


(a . xa . ) 
i D 


matrix given by 





A A 

for i 7 ^ j and with B = (b^j) an (nxm) matrix given by 



(14) 


(15) 


(16) 


Since B as given by Eq. (16) has zero rows except for the 
d^-th, d 2 ~th, . . . ,d m th rows, we need only calculate the correspond- 
ing columns of (si - A - BF) 1 in order to obtain the transfer 

_ 1 /s 1 ^ 

matrix T F (s) = C (si - A - BF) X B = C(sl - A - BF) X B. Moreover, 

BF has zero rows except for the d,-th, d„-th,...,d -th rows and 

^ -a ± z m ^ 

so A + BF is again a block matrix of exactly the same form as A. 
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In other words, A + BF = i- s a matrix of the form 


$ = A + BF = 


$ 

-11 

$ 

-21 


-ml 


-lm 

$ 

-2m 


-mm 


where is a(c^xc^) companion matrix given by 


-li 


0 

0 


1 

0 


d. ,d. n +l ^d. ,d. +2 

l l-l l l-l 


0 

0 


0 

0 


d . ,d . -1 T d.,d. 
i i l i 


and . is a (o^xa . ) matrix given by 


$. . = 
"ID 


0 0 

d) (|) 

, d.,d. ,+l d.,d. 1 +2 

L i' J-l i' ]-l 


d. ,d- , 
i' DJ 


for i / j. These two simple observations are basic to the 
structure theorem (Tl) : 


Tl: Suppose that the pair (A,B) is controllable and 

let T F (s) = C(sl - A - BF) -1 B be the transfer matrix 
of the system x = (A + BF)x + Bw,y = Cx. Then 

( 20 ) 


-1 


T p (s) = CS(s)6 p (s)B m 


( 17 ) 


( 18 ) 


( 19 ) 
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A 

where C 


= CQ 1 * * * f S (s) is the (nxm)matrix given by 


. . 0 

. . 0 


S (s) = 0 


a 2 _1 

0 s Z 0 


a -1 
m 


6 (s) is the (mxm) matrix (6 . . (s) ) with entries given 

— F r f X J — 


S F,ii (s) ■ det<s h. - hi* ZBi { p> ij< s) = 


a ± - 1 

-<f) d.,d. , +l -s< ^d . ,d . -+2 ~ S ^d . , d . for i ? j, and 

i j-1 l j-1 x 3 ' J ' 


B is the (mxm) matrix given by 
-m _ » 


1 b 

c 

0 1 


B m = . . 

“"in ■ • 


0 0 


where B = QB = (b^.) • 
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Proof: In view of the proposition (PI) , we need only show 

A -- - /\ A A -1 A A "I A 

that C (si - A - BF) B = CS (s) <5 p (s)B m . To do this, it will be 
sufficient to show that 


. (si - A - BF) _1 B = S(s)6 1 (s)B_ 

— — — — — “ ■— p in 


(23) 


or, equivalently, that 


A A A 


-1 


(Si - A - BF)S(s) = BB^ 5^ (s) 


(24) 


But Eq. (24) is an immediate consequence of the definitions of 
S(s) and 6 (s) . Thus the theorem is established. 

* — r 

This seemingly innocuous and esaily proved theorem has, as 
we shall see, a number of significant consequences. For a begin- 
ing, we have the corollary (Cl) : 

Cl: Let A (s) = det (si - A - BF) , and C*(s) = CS(s). 

Then, A (s) = det(5 (s)) and if p = m, then 

p ~ p — ■■ 

det T p (s) = (det C*(s))/A p (s) (25) 


where T_(s) = (s) /A_ (s) (i.e. N„(s) is the numera- 

r ^ ~ r — 

tor of the transfer matrix) . 

Proof: By the definition of T„(s), we have T^(s) = 

— r —p 

N p (s)/A p (s). It follows from the theorem that 


N p (s) C*(s)D p (s)B m 


A p (s) det ( 6 p (s) ) 


(26) 


where 6 "'"(s) = D„ (s) /det (6„ (s) ) . However, A (s) and det (6_(s)) 

— r — r “p p — r 

are both monic polynomials of degree (n) and - the entries in N (s) 
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are polynomials of, at most, degree (n-1) . It follows that 
A F (s) = det (6 f (s)) and, hence, that Eq. (25) holds (since 
det(6 T , 1 (s)) = I/det (<$„ (s) ) and det B =1). 

— r "’ll -m 

A A A /N /S >\ 

Now let (A + BF) = A + B F be defined as the (mxn) matrix 

m m m ^ ^ ^ 

consisting of the m-ordered d^-th (k = l,2,...m) rows of (A + BF) , 
and define diag[s CT i] as the (mxm) diagonal matrix whose ii entries 
are s G i. We then have the following corollary (C2) : 

OT • ✓S xs. 

C2: 6 F (s) = diag [ s 1 ] - A m S(s) - § m FS(s) 

This corollary can easily be established by direct substitution 
by using the definitions given above. 

Proof : It then follows that 

a i 

6 (s) = diagfs ] - A S(s) 

-o — m— 

and consequently that 

6 p (s) = <$ o (s) - B m FS(s) 

We note, in particular, that these latter three relationships 

/N /\ 

hinge on the fact that only the m d^-th rows of A and B contain 
the most pertinent information regarding the dynamics of the 
system. The last equation (Eq. (28)) can be established directly 
by using Eq. ( 24 ) . 

We observe that entirely analogous results can be obtained 
for observable systems by a consideration of the dual system 
(refs. 3 and 10) 

x = A ' x + C'v , y ~ § ' x (29) 

which is controllable if and only if Eq. (1) is observable. While 
we will not derive the results for observable systems here, we 
will use them without further ado in the sequel. 


(27) 


(28) 
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2.3 A General Structure Theorem 


Consider the system, Eq. (1), and again let K = [B,AB, . . . 
A n-1 B] . However, we no longer assume that Eq. (1) is controllable 
and so, the (nxnm) matrix K has rank r with r < n. To obtain a 
structure theorem in this general context, we shall consider a 
controllable extension of Eq. (1) and apply the theorem (Tl) . 

With this in mind, we let q = n-r and W be the r-dimensional sub- 
space of R spanned by the columns of K. Denoting the orthogonal 

n j, i 

complement of W by W so that R = W © W and letting 8-, , ... *3 
be a basis of W x , we consider the system 

x = Ax + B e y , y = Cx (30) 

where B is the (n x (m+q) ) matrix given by B = [B,§, ,8 ] • The 

0 c X 

system, Eq. (7) , is controllable and there is a Lyapunov trans- 
formation Q e which carries Eq. (7) into block companion form. We 
note that is a nonsingular (nxn) matrix. It follows that the 
system 

z = Az + Bu , y = Cz (31) 

where A = Q 1 AQ , B = Q B, and C = CQ 1 is equivalent to Eq. (1) . 
0 s ^ e 0 

Moreover, the matrix A is in block companion form, the last n-r 
rows of B are 0, and the lower left-hand block (n-rxr) of A is 0. 
Thus, th© last n-r rows of A cannot be altered by state variable 

/N "* 

feedback of the form u = Fz + w. We now have the following 
theorem (T2) : 

T2 : Let T (s) = C (si - A - BF) _1 B be the transfer 
matrix of the system x = (A + BF)x + Bw, y = Cx. 
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Then 


-1 


T p (s) = 


CS(s)4^ u (s)6 E ; c ( S )B m 

Vu (s) 


."I 


where C = CQ , S(s) is the (nxm) matrix given by 


S(s) = 


1 

s 




0 


0 

0 


0 

1 


0 

0 


0 

0 


’ a 2 _1 

0 s ^ 


' c -1 
m 


, cr-i -1 a -1 

(with b, , Ab-. , ... , A b.. , . . . , A m b a 

- -1 — 1 - -1 - -m — 

"lexicographic" basis of the range of K so that 


m 


i=l 


0 i = r) ' A F,U (S) = dSt -F,u (s) ' -F (s } 


is the (m+q) x (m+q) matrix (6^ ^(s)) with entries given 

by Sp^s) = det (si - i ±i ) .and'« i (s) = -* d , 

- - .01-1. . . ....... - 1 3-1 


s ^d-L^d-j for i ^ j where 


\ - £ "i . »! - 1 

1=1 


( 32 ) 


(3 3 ) 
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for i = m+1 , . 

so that 


A A A 


. , m+q , and A + BF = (cf^.) = [<^-3, 


« F («) = 


6 F,11 (S) •** 6 F,lm (s) 6 F,l,m+l (s) *** 6 F,l,m+q (s) 


6 F,ml (s > 


6_ ( s) 

F ,mm 


6 F,in+. , m+1 ^ 


^F ,m+l ,m+q ^ 


<5^ (s) 1 6 (s) 

-F,c I -F,cu 

l 

i 


^F,m+q,m+l '** ^F, m+q, m+q 


(34) 


e ! %u (s » 

/N 

and where B m is the matrix (mxm) consisting of the 
nonzero row s of B 

Proof : Clearly, we now need only show that 


A A A 1 A 


C(sl - A - BF) 1 B = 


CS ( s ) A (s)<5 l(s)B 
F , u -F -m 


A F,u< s > 


where F = FQ ^ . We shall do this by considering the completely 


controllable system 


• /\ 


Z = Az + B v , y Cz 


(35) 


+6„ (s) involves only constant terms and the off-diagonal terms 

-F,cu 2 

in (s) are constant. 

— F , u 
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with § e = Q e § e and applying the theorem (Tl) . 


/s -1 

Let F = F Q where F = 
-e -e-e -e 

§ e = [§§!'•••»§„]# we have / hy 


so that F - 
-e 

definition of 6 , 
s e 


0 - 


Since 



r i s l 

XV 

B = 

-e 

i 

i 

I 

_ j 

< ffll 


| I 


i q J 


(36) 


and B F = BF. It follows that (si - A - BF) = (si - A - B F ) 

-e-e — - - — - — -e-e 

and, hence , that the transfer matrix of Eq. (35) under the feed- 
back v = F g x + w is given by C(sl - A - BF) However, the 

system represented by Eq. (35) is controllable and thus by 
theorem (Tl) , 


XV XV /\ X\ /V x\ I XV 

C (SI - A - BP) B e = CS e (s)« F (s)B eim+q 


(37) 


where § e (s) is given by 


S(s) | 

0 

1 

1 

I 

1 

-q 


S e (s) = 


and B g is the (m+q x m+q) matrix given by 

B ' 0 

-m | 


B 

-e,m+q 


I 


I 

-q 


(38) 


(39) 


By equating the appropriate blocks in Eq. (37) and noting that 


£;hs) 


(det £F,a (s),aaj £ F,c (s) j ' (adj £f,c (s)) £f , cu (s) (ad 3 £f,u ( s) 
0 f (det 6 F ^ c (s))adj 6 p , u (s) 

I 


det -F,u (s ' 1 det -F,c (s) 


(40) 
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where adj (•) denotes the adjoint of a martix, we deduce Eq. (32). 
Thus, the theorem is established. 


Corollary C3: u (s) is independent of F and the uncontrol - 

lable poles of the system x = (A + BF)x + Bw , y = Cx are 
the zeros of A p u (s) [= Aq, u (s)] 


C3: is simply a statement of the fact that the uncontrollable 
poles cannot be altered by state variable feedback. We also note 
that the factorization (Eq.C32)) involves the well-known pole-zero 
cancellation of the uncontrollable portion of the system (ref. 11). 
Then now state corollaries (C4, C5, and C6) : 


C4: The matrices C , S(s) and B m are invariant und er 
state variable feedback. 

C5 : Let p = m and C*(s) = CS(s). Then the inverse 
system exists if and only if C* (s) is non-singu - 
lar . 


C6 : Let p = m and let Ap(s) = det S F (s). Then 
det(T F (s)) = (det C*(s)T(A F u (s))/A~(s) where 

A F(s) = A F,u (s)A F,c (s) . 

We again observe that entirely analogous results can be ob- 
tained for systems which are not observable by a consideration of 
the dual system (Eq. (29)). We use these results without further 
ado in the sequel . 


3. DIRECT APPLICATIONS OF THE STRUCTURE THEOREM 

Now that the structure theorem has been established, we can 
investigate its application to various questions concerning the 
analysis and synthesis of linear multivariable systems. In parti- 
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cular, in this section we will discuss three "direct" applications 
of the structure theorem. We use the term "direct" to imply that 
the applications are relatively straightforward and require only 
minor modification of the relationships established in Section 2. 

In Section 3.1, we consider the question of realization. In 
Section 3.2, certain questions related to pole assignment via 
linear state variable feedback are discussed. The final section 
(3.3) deals with linear output feedback for pole assignment. An 
example is given in the last section to clarify certain of the 
statements made. 


3 . 1 The Problem of Realization 

The first application we will discuss involves the use of 
the structure theorem to obtain an algorithm for solving the 
problem of realization (refs. 3 and 12). More precisely, we con- 
sider the following . 


REALIZATION PROBLEM: 


Let T ( s ) be a (pxm) matrix whose 

Suppose that 


entries, T. . (s) , are rational functions of s. 

lj 

T ij ^ = n ij ( s ) where n^j (s) and d^ (s) are relatively 

prime and degree n^ (s) < degree d^ (s) . Then determine a triple 
{A,B,C} of matrices such that 


T (s) = C (si - A) 1 B (41) 

and {A , B } is controllable and {A,C} is observable. Such a triple 
is called a minimal realization of T(s) (refs. 3 and 12). 

Ho and Kalman (ref. 12) proved that the realization problem 
has a solution and provided a constructive procedure for determin- 
ing a minimal realization. Here, we present an alternate con- 
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structive algorithm for determining minimal realizations analogous 
to a recent result of Mayne (ref. 13) . A computer program has 
been developed for applying the algorithm. 


The basic steps in the algorithm are: 


STEP 1: Calculation of the least common multiple of 
the denominator polynomials {d. . (s) , . . .d . (s) } in 
each column of T(s) . ^ ^ 


STEP 2: Construction of a standard controllable 

realization {A ,B , C } (not necessarily minimal) . 
c o c 

STEP 3: Construction of a minimal realization by 

applying a suitable transformation to {A' 

c c c 

We shall examine each of these steps, in detail, paying particular 
attention to step 2. 


Now let gj (s) be the least common multiple of the denominator 

polynomials {d.. . (s) , . . . ,d . (s) } (which are assumed, for conven- 

J PI 


ience , to be monic) . 

T* (s) be the (pxm) matrix given by 


Let hj denote the degree of g^ (s) and let 


T*(s) = 


n ll /^i 


n pl (s) 


n !m (s)/g ™ (s) 


m 


n* (s)/g (s) 

pm ' ' 


(42) 


where nt ^ (s) = n^ ^ (s) g^ (s)/d^ ^ (s) . In other words, T*(s) is ob- 
tained from T ( s ) by multiplying each numerator n_„ (s) by 
gj(s)/d^j(s) and replacing each denominator d^^ (s) by g^ (s) . The 
construction of T*(s) completes step 1. 


Since g^ (s) is the least 


m k 

Let n. = £ h . and p. = £ h . . 

1 j=l 3 ** 1 j 

common multiple of {d, . (s) , . . . ,d . (s) } and degree n. . (s) < degree 

-L J P J “*■ J 

d. . (s) and the d. . (s) are assumed monic, we have 
ij 


g. (s) 


h. h.-l 

S-’+Y.tS - 1 +...+Y., 

Jl Dh 


(43) 
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h.-l 

n* . (s) = v, ,.s ^ 

ID i;jl 


h . -2 


+ v . . „s 
132 


+ • • • + V . . , 

i:h. 


(44) 


for all i,j and suitable constants y , v. .. 

J '3k 13k 

companion matrix corresponding to (s) so that 


Let A . be a 
-c,3 


A 

-c r 3 


0 

0 


1 

0 


0 

0 


0 

-Y 


jh_ 


0 

-Y 


jh.-l **• 


1 

“Y 


jl 


and let A c be the (n^xn^) block diagonal matrix given by 

0 


A = 

-c 


A i 
— C,1 


0 


-c , 2 


-c ,m 


(45) 


(46) 


If B is the (n.xm) matrix with zero entries in all but the p.-th 
rows, each of which is zero except for a 1 in the k-th column, 
then the pair (A c ,B c ) is controllable. We now have the proposi- 
tion (P 2 ) : 


C = 

-c 


P2: Let 

C be the 

—0 

(mxn 1 ) matrix given 

by 


V llh 1 

v ll hl -l * 

‘ * V lll 

| V 12h 2 **• 

V 121 | * * * 

v lml 

V 21h x 

• 

V 21h 1 -1 * 

** V 211 

• 

I V 22h 2 *•* 

V 221 ! 
1 

1 

I 

v 2ml 

V TV. 
Plh 1 

V plh 1 -1 

• 

V PH 

1 

1 V p2h 2 

1 

1 

v ! 

p21 | 

v pml 


( 47 ) 


then (A ,B ,C } is controllable realization of T(s) 

• « — C "" C C - — - . . . — , — » — - > ■ — ” 
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Proof : Since {A , B } is controllable, it follows from the 

^ c 

structure theorem (Tl) and the definitions of A ,B ,C that 

-c -c -c 

£ 0 <=I ' = c'lslCwfc,.. (48) 

where B = I , S'hs) = diag [1/g, (s ) , ,1/g (s)] , and C*(s) = 

^ l ^ ~'iii ” c x ni x* 

(n. . (s) ) . Since n* . (s) /g . (s) = n. . (s) /d. . (s) , we deduce that 
-*0 _ n i D D ID 

C(sl-A) B =(n..(s)/d..(s))=T(s). Thus , the proposition 

“G C C 13 

is established. This proposition completes the descriptiorf of 
step 2. 

As regards step 3, we consider the triple {A^C^/B^} and 

apply a similarity transformation Q , of the type used in Section 

“ e n-i-1 

2 . 3, to it. Letting n be the rank of [C^A^C^ . . . A^ C^] and 
setting A^ = Q^Q" 1 , = Q e C^ , B^ = B^q" 1 , we have 

l 
I 

I * 

\ 

I * 

I 
I 
I 

and B 1 = [B'l* ] where C' is (nxp) , A' is (nxn) , and B 1 is 

-c - 1 m,n]_-n - ™ 

(mxn) . Since T(s) = C (si - A ) 1 B , it follows that T’ (s) = 

A A 1a"" ^ OO 

B' (si - A'J C 1 = B' (si - A ' ) ~lc 1 or, equivalently, that 
T ( s ) = C(sl - A X B . Thus, { A , B , C } is a realization of T(s). But 
{A,B,c} is both controllable and observable and, hence, is a 
minimal realization (ref. 12). The triple {A,B,C} is in "observ- 
able canonical form. " The actual available program also produces 
a minimal realization in "controllable canonical form" as well as 
all the relevant similarity transformations. 

3.2 Pole Assig nment via Linear State Variable Feedback 

Several recent papers have dealt with the utilization of 
state variable feedback for closed loop pole assignment in linear 
multivariable systems (refs. 14, 15, and 16). The main result 
established in all of these papers was first presented in 
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reference 16 and is summarized below. 


Let A = ' * ’ * '^n^ ^* e an arl3; *- trar y set of n complex 

numbers A^ , which appear as conjugate pairs whenever Im(A.)^0. 
From theorem (T3) , 

T3 : The pair {A,B} is controllable if and only if, 
for every choice of the set A , there is a matrix 
F such that A + BF has A for its set of eigenvalues. 

The sufficiency portion of the proof was established in 
reference 16 by using a result of Langenhop (ref. 17 ) , However, 
a simple constructive procedure for establishing necessity did 
not appear until Heymann (ref. 15) . The purpose of this section 
will be to present an alternative constructive procedure for 
arbitrarily assigning all n eigenvalues of (A + BF) based on the 
structure theorem. 

The constructive procedure we will employ involves two 
relatively simple observations: (1) the eigenvalues of similar 

matrices are identical (ref. 8) , and (2) the m d, -th rows of 

A A A ^ 

(A + BF) can be completely and arbitrarily specified via F. By 

/A 

(1) , we know that if an F can be found such that all n eigenvalues 

A. A, 

of (A + BF) belong to the set A , then we are finished; i.e. by 

A 

(1) , all n eigenvalues of (A + BF) , where F = FQ, will also belong 

A 

to A. However, by (2) , we note that F can be chosen such that 

A A A 

(A + BF) is a companion matrix with any arbitrary final (n-th) 

A. A. 

row. In particular, we can choose the last row of (A + BF) to 
correspond to the coefficient of the polynominal whose zeros are 
the numbers of A ; i.e., if 

n _ i 

II (s - A.) = s n + a . s n + ... + a. s + a (50) 

. . i n-1 1 o 
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then we can select F so that 


0 1 0 . . . 
0 0 10 


0 

0 


A + BF = 


( 51 ) 


1 


-a 


o 


-a 


1 





A A. 


because all nm terms comprising the m d^-th rows of A + BF can be 
completely and arbitrarily selected by using linear state variable 


feedback. Clearly, the eigenvalues of A + BF are the members 

A 

A and all that is required is an explicit expression for F. 


we define A as the 

rows of A, and A* an 
^-m ^ 

d^-th rows of A + BF 


(mxn) matrix consisting of the m-ordered. 

(mxn) matrix consisting of the m-ordered 
, we observe that 


of 

If 

d. -th 
k 


A 



A A 

+ B F 


-m- 





(52) 


or that 

A -I A. 

£ - £m <53 > 

Necessity in theorem T3 is thus constructively established. 

We note at this point that, in general, Eq. (53) represents 

A 

only one of an infinite number of choices for F , and consequently 
F , for pole assignment in multi-input systems. By altering or 
linearly combining the column vectors of B before transforming 
the system to canonical form via Q, we can achieve any number of 
appropriate feedback matrices. While this might appear to present 
the system designer with an impossible task of selecting the 
"best" F for feedback compensation, we note that this is not 
generally the case in practice for a number of reasons. In gen- 
eral, practical system constraints are present which prevent 
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arbitrary selection of the first feedback matrix F which produces 
a desirable closed loop pole configuration. For example, two 
different feedback matrices may yield the same pole configuration, 
but one of these may be far superior from the point of view of 
reducing system sensitivity to parameter variations. Along these 
same lines, the feedback design which reduces sensitivity might 
be optimal in the sense of a quadratic performance index, while 
the other need not be, although both designs yield the same closed 
loop poles. Constraints might also be placed on the allowable 
magnitudes of feedback gains in order to conserve system power or 
avoid noise amplification. A certain form for the closed loop 
transfer matrix, such as being diagonal (decoupled) , might be 
desirable; in which case, the class of allowable feedback matrices 
would be significantly reduced. If the entire state of the system 
were not directly measurable, one might want to investigate the 
feasibility of employing output rather than state feedback to 
achieve a satisfactory closed loop design. Other constraints, or 
a combination of constraints, can significantly reduce the allow- 
able class of feedback matrices, perhaps to the point where linear 
feedback alone would be unsatisfactory. Most of these points will 
be discussed in more detail in the remainder of this report. 

3 . 3 Pole Assignment via Output Feedback 

At the conclusion of the previous section, we noted that the 
entire state of a given system might not be directly measurable. 

In particular, in many practical systems, various sensors must be 
used to obtain the physical "state" of the system or some portion 
of it. Feedback designs which involve the entire state of the 
system are often difficult to achieve in practice unless addi- 
tional sensors or state estimators (filters) are employed. This, 
of course, involves additional costs which may be avoided if al- 
ternative designs are used. The purpose of this section is to 
discuss a procedure for specifying closed loop pole locations by 
using linear output feedback. The technique we will employ in- 
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volves a direct application of the structure theorem and, in 
particular, corollaries Cl and C2. It should be noted that the 
question we are considering here has a rather elegant solution in 
the scalar case, namely the familiar root locus. Unfortunately, 
there does not appear to be a simple pictorial solution to the 
question of output feedback compensation in the multivariable 
case due to the dimensionality of the problem; i.e. instead of 
feedback from one output variable to one input terminal, we con- 
sider feedback from p outputs to m inputs and, consequently, a 
total of pm individual gains, instead of one. We remark here 
that recent results pertaining to this question have been obtained 
by using a pictorial nested root locus approach (ref. 18). The 
approach taken here, however, will not involve any graphs or 
plots. In particular, consider the multivariable system, Eq. (1) , 
and suppose we wish to calculate the effect on the closed loop 
eigenvalue if output feedback is used, i.e., if 

u = HCx + Gw (54) 

The choice of G does not affect the eigenvalue of (A + BHC) and 
can be neglected. Also, as in the case of state feedback, we will 
consider the canonical equations of motion (Eq. (7) ) rather than 

A A A 

the actual system; i.e., the eigenvalue of (A + BHC). Note that 
H, the (mxp) output feedback matrix is unaffected by a transforma- 
tion of coordinates. The eigenvalue of (A + BHC) and, consequent- 
ly, of (A + BHC) can therefore be determined directly in terms of 
the actual H which is employed. According to corollaries Cl and 
C2, these eigenvalues are equal to the zeros of det (6^ (s) ) , where 

-HC (s) = £ 0 (s) " imSS§ (s) (55) 

The evaluation of the determinant of Eq. (55) directly produces 
the characteristic polynomial of the output feedback state matrix 
(A + BHC) in terms of the pm elements comprising H. 
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This point will now be demonstrated by example. In particu- 
lar, consider the following triple {A,B,C} representing the system 
(Eg. (1)). 



If one transforms the above triple via the similarity transforma- 
tion Q , where 



the resulting triple {A,B,C} is in companion canonical form as 
desired; i.e. 
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and 


C 


0 

0 


0 

0 


1 0 
0 0 


0 0 
1 0 


We note that this system is not completely controllable; i.e. 
A p ^ u (s) = A hc , u (s) = (s + 1) . Also, 


V s) “ 


s J - 2s + 1, 2 


t 2 

-3s , 


s + s + 4 


B 

-m 


1 2 


0 1 


CS(s) = 


s 2 0 


0 s 


We can simplify the computations required to compute the determin- 
ant of 6 RC (s) = 6 (s) - B^HCS ( s ) still further by defining 


i.e. 




2 2 

s + h 1 s - 2s + 1, h^s + 2 

(h 3 - 3) s 2 , s 2 + (h 4 + 1) s + 4 


Solving for the determinant of 6„„(s) in terms of the unspecified 
elements of H* , we obtain 

det 5 RC (s) = s 5 + (h x + h 4 + l)s 4 

+ (h 1 + h x h 4 + h 2 (3 - h 3 ) + 2 ) s 3 

+ (4h 1 - 2h 3 - 2h 4 + 5)s 2 + (h 4 - 7)s + 4 


By inspection, we note that h 4 , h^, and h 3 can be chosen succes- 
sively, in order, to arbitrarily specify the coefficients of s, 

4 2 

s , and s . Furthermore, if h^ ^ 3, h 2 can then be chosen to 
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specify the coefficient of the s 3 term. We also note that the 
constant term is 4, and it cannot be altered by H* ; i.e., the 
product of the five controllable poles of the output feedback 
system must equal -4. Let us carry the analysis to its conclusion 
by assuming that we wish to choose closed loop poles (eigenvalues 
of (A + BHC) ) at s = -1, -2, -1/2, and -3. The fifth controllable 
pole must be -4/3 in order to satisfy the requirement that the 
product of the five controllable poles be -4. The desired closed 
loop characteristic polynomial is thus s 3 + s^ + y^-s 3 + ^p-s 3 + 

p-s + 4. Solving for h^ , h^, and h^ in succession, we obtain: 
h. = 1^- , h, = - — , and h_ = - 599 . Since h- / 3, we can 

4.3± ^68^12 J 3 

solve for h„ to choose rr— as the coefficient of the s term; i.e. 
6094 ^ J 

h_ = ? - ■=- . The actual output feedback gain matrix, H, can now 

^ D O D ^ ^ 

be determined since B is nonsingular; i.e. H = B H* , and in 

-m ^ - -m - 

terms of the elements of H* which we have computed, 

-81.33 41.07 

49.92 -25.33 

This choice for the output feedback gain matrix will yield closed 
loop poles at s = -1, -2, -1/2, -3, and -4/3, in addition to the 
uncontrollable pole at s = -1. 

4. LINEAR OPTIMAL CONTROL 

Over the past few years, optimal control theory has proven 
to be a useful tool for designing linear feedback control systems. 
In this section, we will discuss many of the implications of a 
feedback design based on linear optimal control. We state here 
that few new results will be presented. What we will emphasize 
are certain computational simplifications which can be made over 
previous work, if one formulates this problem in the frequency 
domain and then employs the structure theorem. 

In Section 4.1, we state the optimal control problem consid- 
ered. We then present a known time domain solution and discuss 
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certain characteristics of an optimal feedback design. In the 
next section (4.2) we employ the time domain solution of Section 
(4.1) to derive a frequency domain solution based on the structure 
theorem. We then employ this solution to establish certain fre- 
quency domain relationships implied by an optimal design. In the 
final section (4.3) we discuss a frequency domain solution to the 
optimal control problem by using spectral factorization, and con- 
clude with a discussion of the root-square locus. 

4.1 Problem F ormulation and Prior Results 

The particular problem which we will consider in this section 
is the so-called output regulator problem (ref. 19). That is, 
given the linear multivariable system, 

x = Ax + Bu , y = Cx (56) 

find a control u* , which minimizes the quadratic performance 
index J, where 

J = l ( 5 T C T Cx + u T Ru)dt (57) 

•'o 

The assumptions are that (1) the pair {A,B} is controllable, (2) 

the pair {A,C} is observable, (3) B is of full rank m<n, and (4) 

T 

R is a positive definite matrix; i.e. R = T T for some nonsmgu- 
lar matrix T. The output regulator problem, thus formulated, is 
a time domain (state space) problem. It was initially presented 
and solved by Kalman (ref. 20) for scalar systems, and later ex- 
tended, in part, to include multivariable systems by Anderson 
(ref. 21), who generalized many of Kalman's original results. It 
should be noted that the determination of u* , the optimal control 
which minimizes J, was not the only factor discussed by Kalman 
and Anderson. In particular, as we will show, a feedback design 
based on linear optimal control theory exhibits certain desirable 
properties which can easily be expressed in terms of a frequency 
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domain representation for the system. For this reason, linear 
optimal control theory provides an important link between the 
classical (frequency domain) and modern (time domain) approaches 
to linear system analysis and design. 

We will now summarize the main contributions of Kalman and 
Anderson relative to the stated output regulator problem, and 
later offer an alternative approach and solution via the structure 
theorem. In particular, we state, without proof, the following 
theorem (T4) and two corollaries (Cl and C8) . 

T4 (ref. 20) : Consider the system x = Ax + Bu , y = 

Cx , where (1) the pair (A,B) is controllable, (2) the 
pair {A,C} is observable, (3) B is of full rank m<n, 
and (4) R is a positive definite matrix. The control 
u* , which minimizes J = jf {x C Cx + u Ru}dt, exists 
and is unique, and can be expressed as a linear 
function of the state x ; i.e. 

u* = -R -1 B T Kx (58) 

where K is the unique, positive definite solution to 
* the matrix Riccati equation, 

T -1 T T 

KA + A K - KBR B K = -C C 

Furthermore, the closed loop (optimal) poles of the 

CJ ip 

system are the zeros of det(sl - A + BR B K) , and 
lie in the half-plane Re(s) < 0 ; i.e. the optimal 
system is asymptotically stable. 

C7 (refs. 20 and 21) : If R = I, then 


(59) 



i.e., optimality implies that 

[•]*[I + B T K(jwI - A) _1 B] - I > 0 (61) 

for a ll real co . 

[•]* is used here to denote the conjugate transpose of the 

T —1 

postmultiplying matrix [I + B K(jcoI - A) B] , and > 0 is a short- 
hand notation for nonnegative definite. 

An additional frequency domain relationship implied by an 
optimal design is the following: 

C8 (ref. 20) : Let A(s) = det(sl - A) and A*(s) = 

— IT ~ 

det(sl - A + BR B K) denote the open and closed 
loop (optimal) transfer matrix denominator poly- 
nomials respectively. Then, 



for all real w . 

Equations (61) and (62) represent two important relation- 
ships which all linear feedback designs satisfy. It should be 
noted that Eq. (62) was established for the scalar case only in 
reference 20 (when m = 1) . 

T —1 

The quantity I + B K(sl - A) B represents the return differ- 
ence in multivariable systems (refs. 20, 21, and 22). Furthermore, 
Eq. (61), an expression involving the return difference, has been 
shown to represent a necessary condition for optimality. We can 
now give a heuristic interpretation of Eq. (61); i.e., optimal 
state variable feedback diminishes the effect of plant parameter 
variations (sensitivity) on the closed loop performance of the 
system. This point is covered in detail in reference 20 and will 
not be dwelled on here. We point out, however, that Eq. (61), 
with strict inequality sign, represents a sufficient condition for 
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optimality of a linear feedback design (ref. 21) . The theorem 
and corollaries presented in this section summarize some important 
aspects of linear optimal control. We can now use theorem T4 to 
derive an alternative (frequency domain) solution to the output 
regulator problem. This solution will then be used to extend 
corollary C7 to include cases when R ^ I and to extend corollary 
C8 to include multivariable as well as scalar systems. 

4.2 Optimal Control via the Structure Theorem 

In order to present a frequency domain solution to the out- 
put regulator problem, we will formulate an "equivalent" optimi- 
zation problem in terms of the canonical state z(t) (Section 2.2). 
A frequency domain solution will then be given for the "equiva- 
lent" problem. This solution will then be related to the original 
output regulator problem, defined in terms of the state x(t). 

In particular, if Eq. (7) is used to define the system dy- 
namics, then 


J °° rp^m^ rp 

(z C Cz + u Ru)dt (63) 

0 

~ -1 

where C = CQ and z = Qx. The optimal control, u* (z) , which 
minimizes (Eq. (63)), subject to the constraint (Eq. (7)) is 
unique (ref. 20). Consequently, it must be related to u* (x) , 
as given by Eq. (58) , since the two optimization problems are 
"equivalent." This is indeed the case, as we will now show. In 
particular, we use theorem T4 , and note that in the case of mini- 
mizing (Eq. (63) ) , 


u* = u* ( z) 



Kz 


(64) 


where K is the unique positive definite solution to the matrix 
Riccati equation 


A A /N |-i-| /N ~| /N |-p A /A. pp /N 

KA + A K - KBR B K = -C C 


(65) 
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It can be shown (ref. 23) that K and K (given by Eq. (59)) are 
related through the transformation matrix Q ; i.e. 

K = Q T KQ (66) 

AfTi rp rp 

Substituting Eq. (66) in Eq. (64), and recalling that B = B Q , 
establishes the desired relationship; i.e. 


— 1 T 1 —1 _rp rp^ _ 1 

u*(x) = -R B Kx = -R B Q QKQQ z = 


1 Am A 

-R B Kz = u* ( z) 


(67) 


For convenience, define F* = R '*"B T K and F* = R "'"B^K. Clearly, if 
u*(z) = -F*z can be determined, then u* (x) =~F*x can easily be 

a 

obtained since F* = F*Q. With this in mind, the main result of 
this section can now be stated and established, as follows in 
theorem T5 . 


T5 : Consider the completely controllable and ob- 

■' >4- - - - A " “ A A 

servable canonical system, z = Az + Bu / Y = 

"" oo T 

and performance index J = jT (z C Cz + u Ru)dt, 
where R is positive definite. Let u* (z) ~ -F*z 
represent the optimal feedback control vector 

* " ' " ■ ' " ‘ a A 

which minimizes J . S_.(s) = 6(s) + B F*S(s) 

— -F* - ~m~ 

and then satisfies the relationship (note: 6(s) 

is shorthand for 6 (s) ) : 

o 


m a rp a ~1 rp a __ rp a 7 rp a m a 

61* (-S)B X RB (s) = 6 1 ( - s ) B 1 RB "o ( s ) + S i (-s)C CS(s) (68) 

-F* -m — m -F* - -m — m - - - — 

Proof : We will use theorem T4 to establish this frequency 

domain characterization for the optimal system. 

Rewrite Eq. (65) in a more convenient form — namely 


A A Arp A Afp A A A ArpA 

-KA - A K = C C - KBR B K 


A 

Add and subtract Ks, 


K (si - A) + (-si - A ) K = C C 


KBR 


I a rp a 

B K 


(69) 


(70) 
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and 


Premultiply and postmultiply Eq. 
(si - A) B, respectively; 


/v m ✓s rn _ *1 

(70) by B (-si - A ) 


/N'T 1 aT — 1 /N /N A. T 1 /\ /\ —1 

B (si - A ) KB + B K(sl - A) B = 


/\ rp /Nm — "I /\ rp ^ /v /\ "I /Nm/N /n __ "1 ^ 

B (-si - A ) [C C - KBR B K] (si - A) B 


Now use Eq. (23) , noting that if F = 0 , (-si - A) B = S(-s) 
-1 ^ 

6 (-s)B , and consequently, that 


[ (~sl 



•N rp /\ IT i 

B (-si - A ) 


-1 


a rp __ rp rp 

B^S (-s)S 1 (-s) 


Using Eqs . (15) and (72), we can rewrite Eq. (71) as 


rp m m ✓v./n /% rp ^ _ "I 

B~6 (-s)S (-s)KB + B KS(s)S (s)B m = 

B^6 (-s)S (-s) [C C - KBR B K]S(s)<S (s)B m 

Premultiply and postmultiply Eq. (73) by 6 T (-s)B m T and B m ^6(s) 
respectively, to obtain 


s T 


/\ /N —"I T 1 /N — rp 'P 

(-s)KBB 6 (s) + 6 (-s)B B KS ( S ) = 

m - - -m - — 


rp /\rp/\ — 1 rp /s 

S (-s) [C C - KBR B K ] S ( s ) 


Noting that 

rp _ rp rp rp ^ __ “I ^—1 

[6 (-s)B m T + S (-s)KBT ][TB 6(s) + T B KS(s)] = 

rp /s rp /\ _ ~| rp /N _ m rp 

6 (-s)B RB 6(s) + 6 1 ( -s ) B B X KS(s) + 

-m --m - - -m - -- 

rp ^ 1 rp s\ /N Irp/N 

S (-s)KBB m 6(s) + S (-s)KBR B KS(s) 


( 71 ) 


(72) 


(73) 


(74) 


(75) 
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substitute Eq. (74) into Eq. (75) to obatin 


Hf 1 ^ — T* T T ^ ^ 1 /N — 1 m m /n 

[6 (-slB^T 1 + S (-s)KBT HTB^Sfs) + T i B i KS(s)] = 


*7* 1 rn Am A 

6 (-s)B X RB 6(s) + S (-s)C CS(s) 


( 76 ) 


Now recall that 


-1 -T-T; 


6 p *(s) = §(s) + B m T T "B'KSU) 


(77) 


or 


?i m 1 5 F *(s) = TbJ-^s) + T _T B T KS (s) . (78) 

Substituting Eq. (78) into Eq. (76) yields Eq. (68) and thus 
establishes the theorem. 

Equation (68) represents a concise frequency domain solution 
to the optimal output regulator problem. As we will show in the 
next section, Eq. (68) can be used to directly solve for the opti- 

^ ^ A 

mal feedback matrix F , and consequently for F, via spectral fac- 
torization. Furthermore, variations in the optimal closed loop 
pole locations, as C and R are varied, can easily be determined by 
using Eq. (68) without first solving the entire optimization pro- 
blem. Before considering these applications of the theorem, how- 
ever, we will present an alternative statement and proof of corol- 
lary C7 , and establish corollary C8 for the general multivariable 
case. In terms of the structure theorem and the results derived 
thus far, corollary C7 can be extended to include the case when 
R ^ I. In particular, we have the corollary (C 9) : 

C9 : An optimal feedback design (output regulator 
problem) implies that the following relationship 
holds for all real w . 


[•]*[! + T T B T KS(j0))6 1 (jw)B m T~ 1 ] - I > 0 


(79) 
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Proof ; We premultiply and postmultiply Eq. (71) by T -T and 
T \ respectively. We then substitute the resulting expression 
into the following relationship: 


„rn Am ^vrn ^ rp ^ m ^ >n "1 /n 1 

(I + T B x (-si - A ) - l KBT ) (I + T B K(sl - A) ' L BT X = 

— r P^ r P ^ r P —1 ✓s /n —I -_rn / N r T 1 /N /n — 1 /n 

I + T B (-si - A ) KBT + T B K(sl - A) BT + 


_ p P/\ r P /N'T 1 — 1 /N /\ — 1 /N rj 1 /N /N _ 1 _ ] 

T B (-si - A ) KBR B K(sl - A) BT 


(80) 


We then employ Eq. (23) , which yields the result 


If 


_rp^.m m m /\ /\ "1 _rn Arp a __ 1 

(I + T B^6 (-s)S (-s)KBT )(I + T B KS(s)J 

"] /s rp rp rp a^P a 1 /\ "1 

(I + T B^6 (-s)S (-s)C CS (s)6 (s)B m T ) 

we then let s = jw , we obtain 


(s) B T 

' -m- 


-1 


(81) 


[■]*(! + 


_m /\m ^ _ I 

T B KS(ju)5 


1 i - 


I = [.]*[CS(jw)6 


-1 


(jw)B m T l ] 


(82) 

Clearly, the right side of Eq. (82) is a nonnegative definite 

matrix. The corollary is thus established. When R = I and 
T 

T = T =1/ Eq. (82) exactly corresponds to the results given in 
references 20 and 21. 


We can now establish corollary C8 for the general multi- 
variable case. In particular, we first write Eq. (77) in a more 
convenient form, namely 


6 p *(s) = l m T _1 (I + T~ T B T KS(s) 6 1 (s)B n T 1 )TB in 1 6(s) (83) 

We then take the determinant of both sides of Eq. (83) , obtaining 


— rp H” 1 ✓s — "I -n — 1 

| I + T B KS ( s ) <5 (s)B m T I 
If we now take the determinant of Eq 


\§. F * (s) M£(s) | (84) 

(81), substituting Eq. (84) 
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( 85 ) 


where appropriate and letting s = jio, the resulting expression 
is 


1 ( -jtu ) | | 6 ^,* I 

— = | I + M(u>) | 

| <S T (-jw) | j 6 ( jw) | 

~ -1 ~ -1 

where M(w) = [*]*CS(jio)5 (ju>)B m T , a nonnegative definite 
matrix. Since the determinant of the sum of I and a nonnegative 
definite matrix is always greater than or equal to 1 (ref. 8), 

Eq. (85) implies Eq. (62) directly, and thus establishes corollary 
C8 for the general multivariable case. 


A final point to be clarified in this section concerns cer- 
tain erroneous statements made in reference 20 regarding strict 
inequality of Eqs. (61) and (62) representing a necessary condi- 
tion for optimality in the scalar case. In particular. Condition 
(II) in reference 20 states that a stable control law may be 
optimal only if the return difference is greater than 1 for all 
real u> . The following example refutes this claim as well as 
other erroneous statements which were based on Condition (II) . 

Consider the scalar system, x = Ax + bu , y = cx, where 

r° r°i 

A = , b = I , and c = [0,2]. Suppose we now wish to find 

L-i oj LiJ 

/ 2 2 
(4x„ + u )dt. 

o ^ 

We note that the pair {A,b} is controllable and that R = 1 is a 

scalar. Theorem T4 can therefore be applied directly. The reader 

can easily verify that K, the solution to the matrix Riccati 

equation, KA + A^K - Kbb T K = -c T c is equal to 21. Hence, 

T ~ 

u* (x) = -b Kx = -2 x 2 . According to corollary C7 and, in particu- 
lar, Eq. (61), optimality implies that [•]*[! + B^K(ju)I - A) "'"B] 

- I > 0. For the example, this condition implies that 
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1 _ 2jo) |ji + 2*0 { . ! > o 


or 


l-oo 

4u> 2 


, 1 2,2 = 

(1-0) ) 


1 - 0 )“ 


> 0 for all real oo 


Clearly, 


4u>“ 


/ 2,2 

(l-oo ) 


is indeed >_ 0 for all real oo . However, we note 


that 


1 + 


4(0“ 


,, 2,2 

(l-oo ) 


= 0 at oo = 0, or the return difference < 1 + 


2s 


s 2 +l 


s= joo 


l-oo ‘ 


is equal to 1 at oo = 0 . Hence, Condition (II) , Eqs. 


(34) and (17) of reference 20 fail to hold in general. 


4.3 Spectral Factorization and the Root-Square Locus 

Perhaps the most significant aspect of theorem T5 is that 
u* (x) can be obtained directly from Eq. (68) via spectral factori- 
zation (ref. 24). This fact had been noted previously by Chang 
(ref. 25), Kalman (ref. 20), and Rynaski and Whitbeck (ref. 26). 
However, their expressions, corresponding to Eq. (68) in this 
report, were not in as compact a form as Eq. (68) and, in the 
multivariable case, included the unaltered term (si - A) . We 

note that the right side of Eq. (68) represents the sum of two 
(mxm) matrices. Furthermore, the entries of each of these two 
matrices are only polynomials in the Laplace operator s (no ratios 
of polynomials appear, as in the term (si - A) ^B) . It would 
therefore appear that a spectral factorization solution to Eq. (68) 
would be easier to achieve than prior spectral factorization 
solutions. Indeed, Rynaski and Whitbeck (ref. 26) dismiss spec- 
tral factorization as a "tedious chore" and, instead, adopt a 
"direct" method for obtaining u* (x) . This "direct" method, later 
refined by Whitbeck (ref. 27) , involves some rather cumbersome 
transfer matrix manipulations. The optimal control u* (s) , is 

also expressed in terms of the initial state x(0) instead of F*x. 
Furthermore, their technique involves an initial spectral factori- 
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zation in order to determine the optimal closed loop poles of the 
system. It appears that their "direct" method is not totally 
amenable to automatic computation and involves more information 
than necessary to obtain u* (x) . Nevertheless, it has been demon- 
strated that their method does produce the optimal feedback matrix, 
F*, and the merits of their technique, relative to other solutions 
to the output regulator problem, remain a debatable subject. 

We observed earlier that Eq. (68) represents a compact fre- 
quency domain solution to the output regulator problem. We will 
now establish this fact by demonstrating that 6 *(s) and, there- 

A A 

fore, F* and F* = F*Q can be determined by a spectral factoriza- 
tion of Eq. (68) . Let us define the right side of Eq. (68) as 
$ ( s ) ; 

m ") m Am a 

$(s) = 5 (~s) B m T TB m 6 (s) + S (-s)C CS(s) (86) 

We will now employ two results due to Youla (ref. 24). 

T 

(I) Consider the matrix $(s). If (1) $ ( — j oo ) = $(jto) and 

(2) $(jco) is nonnegative definite for all real co, then $(s) can 
be written as the product of two matrices, W (-s) and W(s); i.e. 

4(s) = W T ( - s ) W ( s ) (87) 

where W(s) and its enverse are both analytic in the half-plane 
T 

Re s > 0 , and W (-s) and its inverse are both analytic in the 
half-plane Re s > 0 . 

(II) Consider that any two solutions, W^(s) and W 2 (s) , of 
the "spectral factorization" of |(s) which satisfy Eq. (87), are 
orthogonally equivalent (ref. 8); i.e. 

W x (s) = UW 2 (s) (88) 

T 

for some orthogonal matrix U (ref. 18) , where U U = I. $(s) , the 
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right side of Eq. (68) , clearly satisfies (1) and (2) of (I) . 
Furthermore, since |£p*(s)| = 0 only in the half-plane Re s < 0 
(i.e., the closed loop optimal system is asymptotically stable 
(refs. 19 and 20), it follows that 6p*(s) and its inverse are 
both analytic in the half-plane Re s > 0 . The matrix TB^S^* (s ) 
thus qualifies as a solution of the spectral factorization of 
4>(s); i.e., if one determines a solution W(s) of Eq. (87) via 
spectral factorization, then by (II) , 

6 p *(s) = B T T -1 yw(s) (89) 

for some orthogonal matrix U. The problem of solving Eq. (68) 

A A A 

for 6p*(s) and, consequently, F* and F* = F*Q, thus reduces to two 
distinct steps: (1) the spectral factorization of $(s), the right 
side of Eq. (68) ; and (2) the implementation of Eq. (89) once 
W(s) , a solution to (1) , has been found. The first of these 
steps is, by far, the more difficult. The spectral factorization 
problem is important in other related research areas, such as 
filtering (ref. 28) and network synthesis (ref. 29), and has 
been the subject of a number of technical reports over the past 
few years. Several investigators have proposed various solutions 
to the spectral factorization problem (refs. 24, 30, and 39), and 
we will not dwell on their relative merits in this report. It 
should be noted, however, that computer programs which perform 
the factorization have been developed (ref. 31) . 

Once W(s) , a solution to the spectral factorization of <$(s) , 

has been obtained, the task of implementing Eq. (89) to determine 

^ _1 . 

6„*(s) remains. In other words, B T U must be determined m 
— F* -m— - 

order to solve Eq. (89) for 6 w *(s). F* can then be determined by 

employing the relations, 6p*(s) = 6(s) + B m F*S(s), or F*S(s) = 

B~^ ( 6 f * ( s ) - 6(s)), and F* = F*Q. The determination of B^T , 

however, is surprisingly simple, and does not depend on the facts 
^ -1 

that B^ and T are already known, and U is an orthogonal matrix. 
We simply recall that, by corollary C2 and Eq. (27) , 
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<5 p *(s) = 6 q (s) + B m F*S(s), where 5 q (s) = <S(s) = diag[s Gl ] - A m §(s) 
Substituting these expressions for S p *(s) in Eg. (89) and solving 
for W(s) , we obtain 


W(s) = U T TB m 1 diag[s 1 ] - (A m - B m F*)S(s) 


(90) 


We now note that the coeffecients of appearing in each of the 

m i-th columns of W(s) are precisely the entries which appear in 
each of the m i-th columns of U^TB^ . Consequently, U TB^ can 
be determined directly from W(s) by inspection, since is the 

highest power of s which appears in each of the m i-th columns of 

T ~-l ~ -1 

W(s) . Once U TB has been determined, its inverse B T U, can 
m „ -in- -' 

then be used to obtain 6 p *(s) via Eq. (89). F*S(s) is then given 
by 


F*S ( s ) = B m 1 [6 p * (s) - 6 (s) ] (91) 

/S 

F* can then be obtained by inspection and F* = F*Q. An example 
which demonstrates these procedures will now be presented. In 
particular, consider the system x = Ax + Bu , y ' Cx, and per- 
formance index _ /• , T„T^, , T„ N , , , 

J=/(xCCx+u Ru) dt , where 

o 



■ 0 

1 


0 

0 

0 " 


0 


0 


0 

-1 


-1 

2 

-3 


0 


3 

A = 

1 

1 


1 

-2 

3 

B = 

-1 

- 

3 


0 

0 


0 

0 

1 


0 


0 


-2 

3 


2 

0 

0 _ 


o 


1 . 


"2 

0 

0 

0 

0 “ 



'2 

o ' 


C = 

0 

1 

0 

0 

0 


R = 





_0 

0 

0 

0 

1 _ 



.0 

1 _ 



We must now find u*(x) = F*x which minimizes J, using spectral 
factorization. The first step is to find a Q which transforms 
this given system to controllable canonical form as represented 
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by Eqs. (7) and (13) through (16). Using the technique outlined 
in Section 2.2 gives 
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3> (s) 


-2s 6 +44s 4 +8s 2 +102, 6s 5 -25s 4 -16s 3 -20s 2 +18s-28 

-6s 5 -25s 4 +16s 3 -20s 2 -18s-28, 19s 4 +5s 2 +8 

$(s) can then be used as the input to Tuel's program (ref. 31), 
modified to yield £ *(s) directly via spectral factorization of 

” _ ^ —i 

$(s) and multiplication by B^T U; i.e. 

s 3 +5.748s 2 +17.764s+31.018, 1.676S-8. 370 

£p*(s) = 

.417s 2 +6. 203S+9.6, s 2 +. 82s-2 . 461 

F*S(s) is determined next by using Eq. (91); 

1. 497s 2 +5.155s-4. 782, -3 . 78 4 s+1 . 01 3 

F*S(s) = 

1. 417s 2 +4.203s+11.6, .82S-2.461 

and F*can then be determined by inspection since S(s) is known. 

-4.782 5.155 1.497 1.013 -3.784 

11.6 4.203 1.417 -2.461 .82 

and since F* = F*Q 

5.155 -3.285 -4.782 1.013 -3.784 

4.203 13.017 11.6 -2.461 .82 

Therefore, F* , the feedback gain matrix associated with the 

optimal system has been determined via spectral factorization. 

— 1 T 

It should be clear that F* = R B K, where K is the solution to 
Eq. (59). 
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The final topic we will discuss in this section is the "root- 
square locus." In particular, it is desirable, in certain appli- 
cations of linear optimal control theory, to readily determine 
the effect on the closed loop optimal poles caused by changes in 
the elements comprising the performance index weighting matrices 
C and R (or T) . A plot of optimal pole location changes caused 
by weighting matrix changes is known as the root-square locus. 

It was first introduced by Chang (ref. 25) to study optimal scalar 
systems. The root-square locus was refined somewhat by Kalman 
(ref. 20) , and later extended to include multivariable systems by 
Rynaski and Whitbeck (ref. 26). In all cases mentioned, the root- 
square locus was a frequency domain tool, and its implementation 
depended on the ability to formulate a solution to the regulator 
problem in the frequency domain and then to find the characteris- 
tic polynomial representing the closed loop optimal system without 
first determining u*(x) ; i.e. F*. Rynaski and Whitbeck (ref. 26) 

appear to be the only investigators who have sought an expression 
for the root-square locus in the general (multivariable) case. 

In particular, they show (Eq. (7-33) in reference 26) that if 
W(s) = C (si - A) _1 B , then 

det(R + W T (-s)W(s)) = F A(S)A (_ 5) ( 92 ) 

where A F *(s) is the characteristic polynomial associated with the 

optimal system; i.e., note that Ap*(s) = | <$ F * ( s ) | = j si - A + BF* | , 

and A(s) is the open loop characteristic polynomial. We point 

T 

out also that the evaluation of the determinant of R + W (-s)W(s) 
involves some rather cumbersome polynomial manipulations which 
can be significantly reduced if the structure theorem is employed. 
In particular, by simply taking the determinant of both sides of 
Eq. (68), we obtain the following relationship: 

|r| A p * (-s) A p * (s) = det (6 T (-s) B~ T RB~^6 (s) + S T (-s) C T CS (s) ) (93) 
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Clearly, Eq. (93) also represents a means of obtaining the root- 
square locus, since the characteristic polynomial of the optimal 
system, A F *(s), is expressed directly in terms of known quantities 
which include the performance index weighting matrices R and C = 
CQ. The latter expression, moreover, is significantly easier to 
evaluate, as the reader can verify. Equation (93) thus represents 
an alternative expression which can be employed to plot the root- 
square locus in the general multivariable case. We finally note 
that Eq. (93) directly reduces to the scalar root-square locus 
expression given by Kalman (Eq. (45) in reference 20) for the case 
when m = 1. 

We will now present an example to illustrate the application 
of Eq. (93) in determining closed loop optimal pole variations 
corresponding to variations in the weighting matrices C and R. In 
particular, consider the system x = Ax + Bu , y = Cx, and per- 
formance index T /*°? T_,T_ , T_. , - . , 

J=/(xCCx+u Ru)dt, where 

o 

B = 

and 




It has been shown that if r > 0, the optimal control law, u* (x) , 

. —1 T 

which minimizes J, is given by u* (x) = -R B Kx = -F*x, where K 

is the unique positive definite solution to the matrix Riccati 

equation (Eq. (59)). The closed loop poles associated with the 

optimal system are therefore equal to the zero of j si - A + 

““1 T 

§5 § 5l = A* (s) . However, Eq. (93) can be used to determine the 

closed loop optimal poles without first solving for u* (x) ; i.e. 
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F* need not be known. For this example, an expression for the 
optimal poles, in terms of the positive scalar r, will now be 
obtained. 


We first note that the system is already in controllable 
companion form (no transformation of state is required) . There- 
fore, by inspection. 


S(s) = 


1 0 
s 0 

Lo i 


Ms) = 


s 0 


-1 s 


CS(s) = I 


B m = I 
-m - 


If we now use Eq. (93) , we obtain 


rA*(-s)A*(s) = det 


2 n 

s -1 


0 -s 


r 0 


0 1 


s 2 0 


-1 s 


1 0 


0 1 


or 


A*(-s)A*(s) = s 4 - s 6 + ~ (-s 2 + 2) 


This expression can now be used to plot the root-square locus by 
using the conventional root locus (ref. 32). In particular, con- 
sider the following block diagram: 



EQUIVALENT BLOCK 
DIAGRAM • 


The root locus (ref. 32) can be employed to plot the variation of 
the closed loop poles of this system with respect to — . Further- 
more, note that the poles of this closed loop system are equal to 
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the zeros of (s 6 - s 4 + — (s 2 - 2)) which, in turn, equal the zeros 

of A*(-s)A*(s). For the example, the optimal poles and their 

mirror images are thus equal to poles of the system shown in the 

block diagram. An approximate root-locus plot of this system is 

given in Figure 1. The arrows indicate the direction of increas- 

2 

ing r; i.e., as r ->- °° , A*(s) -*■ s (s + 1). Note further, that the 
optimal poles display a Butterworth configuration, as noted pre- 
viously in reference 20. 


jcu 



Figure 1.- Root- square locus plot of the example 

5. TRANSFER MATRIX FORMS AND NONINTERACTION 

In this section, we apply the structure theorem to obtain 
some results pertaining to transfer matrix forms. We provide 
answers to certain questions dealing with the noninteraction of 
various input/output pairs. The motivation for this work should 
be rather obvious to those who have dealt with the analysis and 
design of linear multivariable systems. It is almost always 
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■n 


ii min 


desirable to isolate various segments of a multivariable system 
in order to deal with the smaller and less complex subsystems 
individually. Furthermore, a noninteractive or decoupled system 
is usually easier to control. The control literature of the past 
several years contains a rather significant number of articles 
dealing with noninteractive design and additional results still 
appear quite frequently. 

In this section, we will combine certain prior results with 
the structure theorem in order to obtain answers to a number of 
questions dealing with noninteraction which have, thus far, re- 
mained unresolved. In Section 5.1, we will consider various 
questions pertaining to decoupling multivariable systems through 
the use of linear state variable feedback. In terms of the trans- 
fer matrix of the system, decoupling will be used to denote com- 
plete noninteraction of off-diagonal input/output pairs; i.e. 
the transfer matrix of a decoupled system is strictly diagonal 
and nonsingular. In Section 5.2, we extend the results of Section 
5.1 and discuss decoupling via "input dynamics" as well as state 
variable feedback. Roughly speaking, "input dynamics" involves 
the addition of a dynamical system to the input of the system we 
wish to decouple. 

In Section 5.3, a more general type of noninteraction is 
discussed. Here, we are not necessarily interested in complete 
off-diagonal noninteraction, but rather, in achieving a closed 
loop transfer matrix which displays a significant number of zero 
entries. We will present a simple example here and carry it 
throughout the remainder of this section. In particular, consider 
the system x = Ax + Bu , y = Cx, where 


I 

o 

o 

1 — 1 
o 


1 

o 

o 

o 

o 

o 

ro 

1 


o 
1 — 1 

0 0 1 0 1 

B = 

0 0 

1 

1 1 1 -1 -2 

_ -J 


1 — 

o 
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c = 


[“ “1 

L-2 12 oj 


Clearly, the pair {A,B} is in multi-input companion form (Section 
2.1) and, consequently, the structure theorem can be employed 
directly (without) transformation of coordinates); i.e., 

a l - a 2 “ 2 ' 5 m * 


Sts) - 


"l 

1 ’ 


s+1,1 

s 

0 

CS (s) = 


0 

1 


s-2 , 2 

0 

s 




a . 


det(CS(s)) = s + 4 6(s) = diag[s 1 ] - A m §(s) 


or 


6(s) = 


2 

- - 

o 
1 — 1 

1 


2 „ 

s^ 0 


-3 0 0 0 




s +3, 0 





s 0 



2 




i — i 

o 


, 2 

0 s . 


.11-1 -2. 




— s — 1 , (s+1) 



0 s 








Consequently , 


6 1 ( s ) = 


s 2 +3 


(s+1 ) (s 2 +3) (s+1) 2 


-1 T 


The open loop transfer matrix of this system, C(sl -A) B = 


.-1 


CS(s)<$ (s)B m , is given by 


T (s) = CS (s) 6 1 (s) = 


s +2s+2 1 

r 


(s+1) (s 2 +3) ' (s+1) 2 


s 2 -s 


(s+1) (s 2 +3 ) ' (s+1) 2 
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Note that the open loop transfer matrix of this system displays 

complete interaction of all input/output terms; i.e., no elements, 

t. . (s) , of T (s ) are zero, 
il 

5.1 Decoupling via Linear State Variable Feedback 

We now apply the structure theorem to obtain some results 
related to the problem of decoupling using linear state variable 
feedback. This problem has been examined previously by a number 
of authors (refs. 33 through 36) and a number of relevant ques- 
tions have been resolved. More precisely, consider the following. 

DECOUPLING PROBLEM: Let x = Ax + Bu , y = Cx be an m input , 

m output system . Does there exist a pair of matrices {F,G} such 
that the transfer matrix 


C (si - A - BF) 1 BG = T n (s) (94) 

- - - -- -- 

is diagonal and nonsingular? (In other words, does the state 
variable feedback u = Fx + Gw "decouple" the system?) 

A necessary and sufficient condition for the existence of a 
decoupling pair was first given in reference 34. In particular, 
it has been shown that the system 

x = Ax + , y = Cx (95) 


can be decoupled by using linear state variable feedback if and 
only if B* is nonsingular, where B* is the (mxm) matrix given by 


B* 


f l 

c 2 A B 

If 

c A B 

-m- - 


(96) 


with c^, the i-th row of C, and f^ = min[{j: c^A-’b ^ 0}, n - 1] # 
B* and the f^ can also be characterized in the following way 
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(ref. 35): let T p G ^(s) be the i-th row of the transfer matrix 

- F , G (s ) ; then f i ~ =_min tU S ^ + 1 ^F / G,i (s) ^ °^' n “ 1] 

and B*G = limA(s)T_ (s) where A(s) is a diagonal matrix with 

“ “ g->-oo— - F , G “ 

f .+1 

entries s 1 It can be shown (refs. 34 and 35) that B* and the 

f. are invariant under state variable feedback. 
i 

Here, we shall use the structure theorem to answer the 
following questions: 

1. Assuming that the system represented by Eq. (95) 
can be decoupled, what is the maximum number of 
closed loop poles which can be arbitrarily speci- 
fied while simultaneously decoupling the system? 

2. Assuming that the system represented by Eq . (95) 

can be decoupled, which closed loop poles are invariant 
under decoupling state variable feedback? 

3. How can a decoupling pair which specifies the 
maximum number of closed loop poles be implemented? 

While these questions are to some degree resolved in refer- 
ences 34, 35, and 36, we provide a complete and elementary answer 
to them here. 


Let T (s ) be the transfer matrix of Eq. (95). Then T(s) = 
A u (s) ■ -1 

C*(s) - , -y- c ( s )§ m where C*(s) = CS(s) by the structure 

u' ’ -' 

theorem T2. We recall that C*(s) and A (s) are invariant under 

state variable feedback. Now we let p^ (s) be the greatest common 

divisor of the polynomials which are the entries in the i-th row 

C*(s) of C*(s). We note that P^(s) is invariant under state 

variable feedback. We let r^ be the degree of P^(s) and we use 

the notation 8^ to denote the degree of a polynomial (thus, 

r. = 8 ) . We now have the theorem (T6) : 

1 p i 
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T6 : Suppose that the system (Eq. (95)) can be de- 

coupled. Then (1) the maximum number v of the closed 
loop poles which can be arbitrarily specified while 
decoupling is given by 


m 


v 


= 2>i + h + 1( 


(97) 


i=l 


and (2) the invariant poles under decoupling feedback 
m- 

are the zeros of A u (s) and {det C* (s) }/ _ jj (Pj_ (s) . 

Proof: Let {F,G} be any decoupled pair. Then T (s) = 

-1 - - 

C(sl - A - BF) BG is a diagonal matrix with entries n^^ (s) /d^^ (s) 

where n^^(s) and d^^(s) are relatively prime. We note that, 


-I U !_ < , -] 

since f. = min{j:lim s^ T_ _ . (s) ^ 0},* 3 = 

1 “r # vji . i — n . . 

s+°° - - 11 


= 3, - f. - 1, 

d. . i 

11 11 


It follows from corollary C6 and the definition of the p^(s) that 


m n. . (s) m A (s) 

* 5^= . n P i (s)det C fl (a) ^ 

1=1 li 1=1 F 


A (s) 


det G 


(98) 


where C*(s) is the matrix with rows C*(s) = r—y 

11 11 1 

A p (s) = A u (s)A p c (s ), we have 


C* (s) . Since 


m 


F,c 


■E 


(r. + f. + 1) 

l l 


3* 

II 


i=l 


(99) 


where 3* is the degree of det C* (s) and 3_ is the degree of 

11 —11 r_ / C 

A„ (s). Now, it is clear from theorem T2 that 
F , c 


^F,G,i (s) S' 1 i; 1 «F,c (s) = Si (s > 


( 100 ) 


and, hence, that n^^(s) is a common divisor of the entries in 


*Note that B* is nonsingular. 
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C*(s) (since n^^(s) and d^^(s) are relatively prime). In other 

words, n. . (s) must divide p. (s) and so 3 < r. . Since no more 

11 n . . = i 

m 11 m 

than 53 9 J poles are assignable through {F,G} and 53 = 

i=l a ii i=l ii 

m m 

53 (9 + f . + 1) , we deduce that at most v = 53 ( r - + f. + 1) 

i=l n ii 1 i=l 1 1 

poles are assignable while decoupling. 

Writing T„ _(s) as a diagonal matrix with entries q. . (s)/A„(s) 

~r f Lr Hr 

n ii (s ) /dii (s) , we deduce that q^^(s) must divide p^(s)Ap(s) or, 
equivalently, that 


q ii (s) P ± (s) 

A p (s) q ± (s) 


( 101 ) 


for i 


1 


. . ,m and polynomials q^(s) 

m 


with 

m 



r . + f . + 1 . It 
x 1 


follows that det T n „(s) = II (p.(s))/ II (q.(s)) and, hence, from 

™r / O • -l 1 n 1 

- - 1=1 1=1 

Eq. (98) that 


m 

A (s) = det C* (s)A (s)det G II q. (s) 

F -II u - i=1 x 

det C* ( s) m 

= — — A u (s)det G II q^(s) (102) 

n p. (s) 1=1 

i=l 1 

Since C*-j.(s) is invariant under decoupling feedback, it follows 

that the zeros of A (s) and det C*(s) are invariant poles under 

u -II 

decoupling feedback. 

Thus, to complete the proof we need only construct a de- 
coupling pair {F,G} such that the resulting polynomials q^(s) are 
arbitrary polynomials of degree r^ + f^ + 1. To begin with, we 
note that transfer matrix 
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Illllllllll I lllllllll 


A u (s) -1 - A u (s) -1 

T(s) - C*(s) «o, c (s) Sm = t(s)S! x (») sVr «0,o (s >S 

u - U 


m 


where P(s) is a diagonal matrix with entries p^(s). Setting 


-ii ( s ) - £ll (s) r4?r Cc (s)S m 

u 


we can easily see that r. + f. = min{j:lim s-* + '*'T T . (s) ^ 0} 

1 1 ll/i 


(103) 


and 


S-»-oo 


that B* = limA (s)T (s) = B* where A (s) is a diagonal matrix 

g->-oo 

r.+f .+1 

with entries s 1 1 (note that the p. (s) are monic) . Moreover, 

A 

as C*(s) is given by CS(s) and p^(s) is the greatest common divi- 
sor of the entries in Ct(s) , we can write C*(s) in the form 

/V 1 A J L 

Cj^Cs) for some constant matrix C (where S(s) is given by Eq. 

(33)). In other words, T^^(s) is the transfer matrix of the 
system x = Ax + Bu , y = C^x where = C^Q and Q is the 

Lyapunov transformation corresponding to Eq. (95). Since P(s) is 
diagonal it will be sufficient to construct a decoupling pair 
{F,G} for the system 


x 


Ax + Bu , Y n 




(104) 


such that the closed loop poles are arbitrarily placed. However, 
letting d^ = r^ + f ^ and applying the synthesis procedure of 


reference 34 (p. 655) , we find that Eq. 


and all its closed loop poles assigned. 

d.+l i 

i r> _i_D 


(104) can be decoupled 
d 


To be more explicit, if q^(s) = s 
coupling pair is given by 


" E 11 ^ 3 ' 

j=0 D 


then the de- 


F = B* 1 [ E M,C TT A k - A*] , G = B* 1 (105) 

0 

where d = max d^ , the are diagonal matrices with entries m^ 
*Clearly , it is enough to consider the case of a monic q^(s). 
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) (i.e. the i-th 


(i.e. M k = diag [m*, . . . ,m£] ) and A* = * 


d.+l 


d.+l 


row of A* is given by C TT -A 1 ) . This completes the proof. 

***’ “ -L X f 1 


We are now in a position to determine whether or not the 
system presented at the beginning of this section can be decoupled 
via linear state variable feedback. If the answer is yes, several 
related questions can also be resolved. We first note that B* , 
as defined by Eq. (96), for the example isQ q and, consequently, 
is singular. Hence the system cannot be decoupled by using linear 
state variable feedback alone. This provides some motivation 
for the next section. 


5.2 Decoupling via Input Dynamics 

In the previous section, we considered the question of de- 
coupling, using linear state variable feedback. We may now raise 
some pertinent questions related to decoupling via state variable 
feedback. In particular, suppose B* were singular. Would it 
still be possible to decouple the system by using some alternative 
technique? Also, suppose that certain of the poles which could 
not be altered by decoupling state variable feedback were unstable. 
Could the system still be decoupled and simultaneously stabilized 
by using an alternative decoupling technique? The purpose of this 
section is to provide answers to these two questions. 

Thus far, we have defined decoupling in a somewhat restricted 
sense, namely decoupling via linear state variable feedback. 
Recently, it has been shown that it is possible to decouple a 
larger class of systems than those which can be decoupled via 
state variable feedback alone, through the use of "input dynamics" 
(ref. 37) . Roughly speaking, we will now permit the addition of 
dynamics to the input of the system (replace G in the closed loop 
transfer matrix by G(s), a k-th order dynamical system). It will 
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be shown that certain systems which cannot be decoupled via state 
variable feedback can be decoupled by appropriate selection of 
G(s). In particular, the example presented at the beginning of 
this section falls into this category. Furthermore, in certain 
cases, poles which are invariant under decoupling state variable 
feedback can be altered if an appropriate G(s) is employed. It 
should be stressed that decoupling via input dynamics involves 
state variable feedback as well as G(s) in the most general case. 
This point will be clarified in the remainder of this section. 

We start by defining the class of systems we will consider 
and exactly what we mean by "input dynamics." In particular, we 
will consider systems whose dynamics can be expressed by Eq. (95). 
B and C are assumed to be of full rank m, where m n. 

We next define "input dynamics" as the k-dimensional linear 
system whose dynamical equations can be expressed in the time 
domain as 


g = Kg + Lv 


u e = Mg + Nv 


(106) 


where g is a k vector called the state, v an m vector called the 
input, and u g an m vector called the output. K, L, M, and N are 
constant matrices of the appropriate dimensions. The transfer 
matrix relating the Laplace transform of the output, u (s) , to 
the Laplace transform of the input, v(s) , of Eq. (106) will be 
defined as G(s), the input dynamics; i.e. 


G ( s ) 


= M(sl - K) 


L + N 


(107) 


We can now define a composite system by combining systems (Eqs. 
(95) and (107)) via the relationship 


u = Fx + u 
- — -e 


(108) 
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The composite system of total order (n+k) can now be expressed in 
time domain as 


X 

= Ax 

+ 

Bu 

Y 

= Cx 



u 

= Fx 

+ 

u 

-e 

2 

= 52 

+ 

Lv 

u 

= Mz 

+ 

Nv 


e 


If we define x c = [gj as t ^ ie coin P os i te state, Eq. (109) may be 
written more succinctly as 



- - 


- 



- - 


X 


A+BF, 

BM 


BN 

x = 


— 



x + 


-c 





-c 



2 


o 

K 


L 


Y c = Y = [C,0]x c 


( 110 ) 


Pictorially, we have the following block diagram representation 
(Figure 2) for the composite system. 



u e (s) = (M (si - K) 1 L + N}v(s) y(s) = CS (s) Sp 1 (s) B m u e 
Figure 2.- Block diagram of composite system 
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Clearly, the composite transfer matrix T (s) relating the output, 

“ c 

y(s), to the new input, v(s), can be written as 

~ -1 A u ^ 

T (s) = CS(s)6 X (s)B G(s) A U (s , (111) 

u v ’ 

where G(s) is given by Eq. (107). 

We can now define the concept of "decoupling via input dy- 
namics" as follows. Consider the composite system, Eq. (110). 

We will say that the system (Eq. (95)) can be decoupled via input 
dynamics if there exists a pair of matrices {F,G(s)}, where G(s) 
is given by Eq. (107) , such that the composite system, Eq. (110) , 
can be decoupled via linear state variable feedback. 

Some observations will now be made. First, it should be 
noted that if system, Eq. (95) , can be decoupled by using linear 
state variable feedback (if B* is nonsingular) , then input dynamics 
are not required to decouple the system. However, we recall that 
certain eigenvalues of (A + BF) are invariant under decoupling 
state variable feedback, namely the zeros of A u (s) and det(C*^(s)). 
These eigenvalues (poles) appear in "cancelled" pole-zero pairs 
in the final decoupled transfer matrix. If one or more of these 
poles lie in the half-plane Re s>0, the resulting decoupled system 
will be unstable. Therefore, although B* is nonsingular, decoupl- 
ing via linear state variable feedback may produce an unstable 
system. One might therefore ask if it is possible to employ in- 
put dynamics in these situations in order to decouple the system 
and also ensure stability. The answer to this question is yes, 
provided A u (s) is Hurwitz (i.e., the unstable roots are zeros 
of det (C* (s)) only). Clearly, if A^(s) is non-Hurwitz, any 
linear feedback system will fail to achieve stability. These 
points will be covered in more detail in the remainder of this 
section . 

A pictorial representation of the definition of decoupling 
via input dynamics is now in order. In particular, according to 
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the definition, linear state variable feedback may be employed 

twice in order to achieve decoupling via input dynamics. First, 

we seek a pair {F,G(s)}, which defines the open loop composite 

system (Eq. (110) ) , a system which can then be decoupled via 

linear state variable feedback. Once the pair {F,G(s)} is chosen, 

an additional pair {F ,G } can then be used to actually decouple 

the system (Eq. (110)). The employment of the second pair, 

{F ,G }, is not always required, as we will show. Figure 3 is a 
" c c 

pictorial representation for decoupling via input dynamics by 
modifying Figure 2 to include the additional pair {F ,G }. The 
subscript c denotes composite feedback and feedforward for the 
(n+k) dimensional system (Eq. (110)). 



Figure 3.- Decoupled composite system 

Summarizing, we note that decoupling via input dynamics 
involves two steps: (1) selection of a pair {F,G(s)} which in- 

creases the dimension of the state, resulting in a composite 
system of dimension (n+k) , which can then be decoupled via compo- 
site linear state variable feedback, and (2) selection of the 
constant decoupling pair {F c ,G c } which then decouples the compo- 
site system. We also remark that once step (1) has been accom- 
plished, the techniques outlined in Section 5.1 can be employed 
to accomplish step (2) . 
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We will now characterize the class of systems (Eq. (95)) 
which can be decoupled via input dynamics. Clearly, this class 
must include those systems which can be decoupled by using linear 
state variable feedback. Fortunately, this class can be charac- 
terized readily by using the structure theorem. Furthermore, we 
can also say something about the stability of systems which can be 
decoupled via input dynamics. In particular, we can now state and 
prove the following theorem (T7) : 

T7 : Let x = Ax + Bu , y = Cx denote an m-input, m-output 
system whose transfer matrix is (in terms of the struc- 

X *“"I * 

ture theorem) CS(s)S (s)B A (s)^A (s) . There exists a 

— - -mu u ______ 

pair of matrices {F,G(s)} such that this system can be 
decoupled via input dynamics if and only if CS(s) is 
nonsingular; i.e. if and only if the system is invertible. 
Furthermore, if A u (s) is Hurwitz, the pair (F,G(s)} can 
always be chosen to ensure the existence of a pair {F c ,G c ) 
which decouples the composite system and simultaneously 
ensures that all (n+k) closed loop poles lie in the half- 
plane Re s<0. 

Proof ; In order to establish necessity of the first state- 
ment, we will assume that a suitable pair (F,G(s)} exists. The 
transfer matrix of the open loop composite system is then given 
by Eq . (Ill) ; i.e. 


T c (s) 


CS(s)6 F 1 (s)B m G(s) 


A 

u 


(s) 

W) 


T c (a) represents the transfer matrix of a system which can now 
be decoupled via linear state variable feedback. Consequently, 

T (s) is nonsingular and invertible. Clearly, this implies that 

A ^ 

CS(s) is also nonsingular and invertible. In order to establish 
the fact that the nonsingularity of CS(s) is sufficient to achieve 
an asymptotically stable decoupled system via input dynamics, a 
constructive algorithm will be employed. Since the algorithm is 
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rather involved, no interpretation of the various steps taken will 
be given until the theorem has been established. At that time, 
we will demonstrate its application by example and comment on the 
more important steps which were employed. The steps will be num- 
bered for convenience and later interpretation. 

STEP (1) : We begin by considering T (s) , the transfer matrix 

of the composite system. We then factor CS(s) as in Section 5.1; 

i . e . 


CS (s) = P (sJCJj (s) (112) 

where P(s) is a diagonal matrix consisting of the ordered greatest 

A 

common divisors, p^(s) , of each row of CS(s). We now consider the 
altered transfer matrix, T^ (s) , obtained by factoring P(s) out 
of T c (s) or, equivalently, premultiplying T c (s) by P ^(s); i.e. 


T ( s ) 

-II 


-II (S) -F (s) § m - (s) A U (s) 

u 


(113) 


STEP (2) : We will now obtain a pair (F,G(s)} which diagon- 

alizes Eq . (113). Solving Eq . (113) for B G(s), we obtain 

i m G(s) = « < S )<C* l(S ))- l T < S ) ^4iT <114) 

-II U 

Now, (C*j(s)) = adj (C*-j. (s) t | C*^ (s) | , where ad j ( • ) denotes the 

adjoint, and |*| the determinant of a square matrix. Note that 

/N 

this step is possible since CS(s) was assumed to be nonsingular. 

STEP (3): Next, we factor Ic^fs) | as the product of two 

monic polynomials, p g (s) and p u (s), and some constant c, where 

p (s) and p (s) represent the stable and unstable (or conditionally 
s u. 

stable) zeros of |c|j(s)| respectively; i.e. 

\C* X (s) | = cp s (s)p u (s) (115) 
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Substituting this expression into Eq. (114) , we obtain 


G (s ) 


(s) adj (C*! (s) )T C (s)A u (s) 
' cp s (s) P u (s) A u (s) 


( 116 ) 


STEP (4) : We now set S„(s) equal to a diagonal matrix whose 
a -i - 

entries are (s+A) x ; i.e. 


a . 

6 f (s) = diag [ (s + A) 1 ] (117) 

where A is a positive constant equal to one of the zeros of p (s) , 

s 

if possible. This choice for 6„(s) ensures the asymptotic stabi- 

“ r 

lity of the original portion of the composite system, provided 
A (s) is Hurwitz. Simultaneously, a certain amount of pole-zero 
cancellation will be possible' in reducing Eq. (116) to simplest 
terms. Pole-zero cancellation is desirable if we expect to keep 
k, the dimension of G(s), as low as possible. 

STEP (5): Since G(s) must also represent the transfer matrix 
of a stable system, T£_^(s) WX H now t> e set equal to a diagonal 
matrix which cancels the unstable polynomial, p (s) , now appearing 
in the denominator of G(s); i.e., we let 


Tp (s) 
-II 


diag 


cp u (s) 

^i 

(s+A) 


(118) 


The integers, q^ , i = l,2,...,m, have yet to be determined. 

STEP (6) : Equations (117) and (118) are now substituted into 
Eq. (116) in order to obtain G(s) in terms of the unspecified 
integers, q^ . G(s) can then be written as 

'"-1 

G(s) = B diag 

- -m 


a . 

(s+X) 

adj (C* x (s) ) diag 

cp u (s) 

P„ (s) 

q i 

s 


(s+A cp u (s) 


A (s) 
u 


A (s) 
u 


( 119 ) 
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STEP (7): The ' s are chosen as small as possible, consis- 
tent with the requirement that G(s) represent the transfer matrix 
of a physically realizable system, namely Eq. (106). The actual 
input dynamics we use to decouple the system can now be chosen to 
be a minimal realization of the transfer matrix (Eq. (119)). The 
algorithm has now been carried far enough to establish the theorem. 
In particular, we note that after all pole-zero cancellations are 
made in Eq. (119), the only poles remaining are either equal to 
-A, or are zeros of the Hurwitz polynomials p (s) . Since j 6 F ( s ) | 
and A^(s) are Hurwitz polynomials, the latter by assumption ,” the 
composite system (Eq. (110)) defined by Eq. (117) and a minimal 
realization (Section 3.1) of Eq. (119) is asymptotically stable. 
Furthermore, the transfer matrix, T c (s), of the open loop com- 
posite system is given by 

T c (s) = P(s)diag 


cp u (s) 

(s+A) 


( 120 ) 


a diagonal matrix (a decoupled system) with poles in the half- 
plane Re s<0. The theroem is thus established. 


In general, we do not want a multiple root at s = -A. 

Rather, we would like to arbitrarily specify as many closed loop 
poles of the decoupled composite system as possible. This is 
where the composite state variable feedback, {F c ,G c }, enters into 
consideration. The open loop system, obtained by selecting F and 
G(s) according to the algorithm just outlined, is already decoupl- 
ed. Therefore, composite state variable feedback, {F ,G } can now 

m c c . 

be used to adjust (Eq^) °f the closed loop poles of the composite 

system (see Section 5.1). 


The example presented at the beginning of this section can 
now be used to illustrate this point and the various steps em- 
ployed in outlining the algorithm. In particular, we noted at 
the conclusion of Section 5.1 that the example represented a 
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system which could not be decoupled by using linear state variable 

/N 

feedback alone. However, for the example, we note that |cS(s)| = 
|cS(s)| = (s + 4); i.e., since CS(s) is nonsingular, input dynamics 
can be employed to decouple the system. We will now go through the 
seven steps outlined in the algorithm in order to derive a feedback 
pair, {F , G ( s ) } , which decouples the given system. 

STEP (1): No further factorization of CS(s) is possible; 
i.e. P (s) = I and C*j (s) = CS(s). 

STEP (2); Here, we calculate the inverse of C|j(s), or CS(s) 
in the case of the example under consideration: i.e. 


[CS (s) ] 



s+4 


adj (CS (s) ) 
CS (s) 


Furthermore, note that in the example A (s) =1 and B = I. 

u -m - 

Hence, according to Eq. (114) , 


G (s ) = 6 (s) [CS (s) ] 1 T r (s) 

-F " -Cjj 


STEP (3) : For example, p (s) = (s+4) and p (s) = c = 1. 

S LI 

This step is necessary so that any unstable roots of |c*.j.(s) | 

can be isolated and later cancelled by T_ (s) . 

-II 

STEP (4) : For the example, the choice of X is obvious; i.e., 

since p (s) = (s + 4) , X should be chosen equal to 4 and 
s 

(s+4) 2 , 0 
0 , (s+4) 2 

This particular choice for <$ F (s) accomplishes two things. First, 
it partially assures us of a stable composite system. It also 
results in eventual pole-zero cancellations in the computation of 
G(s) , the added dynamics. This maintains k, the dimension of 


fip(s) = 


66 



G(s), at a relatively low value. 

STEP (5): We now employ Eq. (118); i.e. 


T (s) = 

-II 

If P u (s) ^ 1, this step would have resulted in pole-zero can- 
cellation of unstable roots which otherwise would have appeared 
in the denominator polynomial of G(s). 

STEP (6): We now employ Eq. (119) to obtain a frequency- 
domain description of G(s); i.e. 


G (s ) 
or 


STEP (7) : We now select the q^'s as small as possible, con- 
sistent with the requirement that G(s) represent a physically 
realizable system. By inspection, this implies that q^ = q^ = 2 
or that 
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Note that with G(s) equal to the above asymptotically stable 
system and 


5pfs) 


(s+4) 2 0 

; i.e. CS(s) 6 1 (s)B = 

s+1 

1 

0 (s+4) 2 

F 

s-2 

2 


( s+4 ) 2 


another asymptotically stable system, the open loop transfer 
matrix of the composite system, CS(s)6_ 1 (s)B G(s), is 

“ r “IU"" 


T r (s) = T p (s) = 

-II - 


(s+4) 


(s+4) 


a result which can be obtained either by performing the multi- 
plication indicated above or simply employing Eq. (118). 

We will now continue our analysis of this example by comput- 
ing k, the dimension of a minimal realization of G(s) and then 
describe how composite state variable feedback can be used to 
arbitrarily reassign certain of the closed loop poles of this 
decoupled system. 

Employing the results of Section 3.1, we can readily obtain 
a minimal realization of G(s); i.e. 

G ( s) = M(sl - K) -1 L + N 


where 


a scalar 



K = -4 


L = [2 -1] 
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and 


r o 


N = 


-1 


0 

1 


as noted previously. G(s) thus represents the transfer matrix of 
a first order system; i.e. k = 1. A time domain representation 
for the composite decoupled fifth order system can now be written 
directly : 


and 


A = 
-c 


A+BF , BM 
0 , K 


0 1 | 0 

-16 -8j_ 0 
0 

0 

0 0 0 


0 

0 


I 


0 j 0 

0 ! -16 -81 3 


I 


0 I -4 


0 0 



c c = [C 0] 


1 

-2 


1 1 
1 2 


0 I 0 

I 

I 

0 ! 0 


Using the results of the previous section (5.1), we can 
easily verify that this composite fifth order system can be de- 
coupled via linear state variable feedback; i.e. 


* 



CAB 


-c-c-c 


I 
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Consequently, four of the five closed loop poles of the decoupled 
system can be arbitrarily specified via the algorithm given in 
reference 34 and repeated in Section 5.1 (Eq. (105)). It can also 
be shown that the single pole which cannot be altered while de- 
coupling is s = -4. However, this pole will not appear in the 
final closed loop transfer matrix because of pole-zero cancella- 
tion . 


5.3 Two Noninteraction Algorithms 


In this section, we will outline two algorithms which can be 
used to achieve a degree of noninteraction in linear multivariable 
systems. We have just presented a fairly comprehensive review of 
decoupling (complete off-diagonal noninteraction) . In certain 
systems, however, we might not be willing to add dynamics in order 
to achieve decoupling or to stabilize a system which can be de- 
coupled via linear state variable feedback. However, we might 
want to zero certain closed loop inout-output pairs, if possible, 
using linear state variable feedback. The structure theorem can 
be used to provide partial answers to questions involving the non- 
interaction of various input-output pairs. The reader will note 
that this section lacks the formalism of previous sections. We 
will merely present two design algorithms for noninteraction and 
demonstrate them via the example presented earlier in Section 5. 
Both of the algorithms employ the structure theorem and, conse- 
quently, depend on the ability to express the closed loop trans- 
fer matrix of a. system as 



~ -1 ^u^ s ^ 

£§<s)« f < = >B m G rl¥) - 

u 


( 121 ) 


We still assume that the system has m linearly independent inputs 
and m outputs. 

The purpose of the first algorithm (single row altered) is 
to demonstrate a technique which can be employed to achieve almost 
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diagonal noninteraction (decoupling) . The exception is that a 
certain row of the closed loop transfer matrix will not contain 
just one nonzero entry but could, in general, contain all nonzero 
entries. To demonstrate this algorithm, we observe the effect of 

A 

altering a single row of C; i.e., we let 

C = C + R (122) 

-q 


where R is an (mxn) constant matrix with (m-1) zero rows and only 
one nonzero row (the q-th row) . We would like to choose R in such 

A 

a way that CgS(s) has certain desirable properties from the point 
of view of the decoupling Section 5.1. In particular, we would 

/N 

like to choose R such that (1) the polynomial |c S(s) | is Hurwitz, 

and (2) the matrix B* (defined in the same way as B* (Section 5.1) , 

but with C replacing C) is nonsingular. If this is possible, the 
• ^3 -A 

system x = Ax + Bu , = C^x , where = QC^ , can be decoupled 

via linear state variable feedback; i.e., by appropriate choice of 

the pair (F,G>, C S(s)6 (s)B G will be diagonal and nonsingular, 

— — — rr— —r — m — 


and all zeros of ?6 (s)]" will lie in the half-plane Re s<0. 

“r 

Now, by Eq. (122), the closed loop transfer matrix of the 
actual system we are dealing with can be expressed as 


CS (s) 6 1 (s) B G = C S (s) 6 1 (s) B G - RS (s) 6 1 (s) B G (123) 

— r -m- ""g - “r “ill - ~r_ ill 

— 1 ^ 

where RS(s)6„ (s)B G is an (mxm) matrix of transfer functions 
— -F -m- 

with zero entries everywhere but the q-th row. Consequently, 
the closed loop transfer matrix of the given system displays 
complete off-diagonal noninteraction with the exception of the 
q-th row. Our design objective has thus been met provided R can 
be suitably chosen. 

We will now demonstrate this algorithm, using the example 
given at the beginning of Section 5. In parituclar, in the ex- 
ample , 
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a nonsingular matrix. The altered system can therefore be de- 
coupled via linear state variable feedback. Two of the four 
closed loop poles of the decoupled system can be arbitrarily 
selected (the other two are s = -1 since |C S(s) | = (s+1) ) . If 

we select both of these poles at s = -2 , Eq. (105) can then be 
used to find the required feedback pair {F,G}; for this example, 

F = B* -1 [M C - A*], where 
-q L -o-q 
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0 1 


-2 0 


31/4 -1 1 5/2 _ 

B* 1 = 


, M = 


, A* = 


-q 

-1 9/4 

-o 

0 -2 


-3 -202 


Solving for F, we obtain 


F 


70 -4-2 

8 7/8 -7 -4 


Therefore , 


5 p (s) = <5 (s ) - § m FS(s) = 


s 2 -4 


2s+4 


- 9/2 s-9 s +6s+8 


9 9 — 1 

and | 6 ( s ) | = (s+1) (s+2) as required. Since G = B* , we can 

_ rl 

write the closed loop transfer matrix of the altered system; 

1 


C S (s) 6* (s) B B* 1 
-q- -F -m-q 


s+2 


0 


0 


s+2 


(s+1) 

(s+1 2 


We now use Eq. (123) to determine the closed loop transfer 
matrix of the actual system we are dealing with; 

r 2 2 

s z +s-3 , -s z - 11/4 s + 1/2 


RS(s)6 p 1 (s)B m G = 


0 


0 


J s+2 


(s+1) (s+2) 


s+2 


and 


CS (s) S- 1 ( s ) B B* 1 
— — F -m-q 


[ 


s+4, s 2 + 11/4 s - 
0 , (s+1) 2 


1/2 J 


(s+1) (s+2) 


s+2 

s+2 
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We have thus succeeded in obtaining noninteraction in all 
off-diagonal elements of the closed loop transfer matrix with the 
exception of those in the first row. 

ALGORITHM 2 (Row Combination) 


The second algorithm is a modification of the first. It 

A 

also relies on an alteration of C and utilization of the results 
in Section 5.1. To demonstrate this algorithm, we observe the 

effect of premultiplyiing C by a nonsingular matrix P; i.e., 

/\ 

linearly combining various rows of C. In particular, we let 


C = PC (124) 

P 


where P is a nonsingular (mxm) matrix. We would like to choose 

P in such a way that the system x = Ax + Bu, y = C s, where 
a p p 

C = QC , can be decoupled via linear state variable feedback. 

"P — P 

In other words, we would like the matrix B* (defined as in 

-P 

Section 5.1, but with C replacing C) to be nonsingular. Suppose 

P 

we can decouple the altered system via linear state variable feed- 
back. This means that the closed loop transfer matrix of the 

A # ~j A 

altered system, PCS(s)6 (s)B G, will be diagonal and nonsingular. 
We obtain the closed loop transfer matrix of the actual system 
under consideration by premultiplying PCS(s)6 (s)B G by P 

1 -t ITT 

Note that the nonzero entries of P will produce corresponding 
nonzero entries in the closed loop transfer matrix of the actual 
system, since P premultiplies the diagonal transfer matrix. 

Hence the choice of P, and therefore P ^ , directly determines 
whether the various entries of T_ „(s) will be zero or nonzero. 

-i? , 


We will now demonstrate this algorithm, using the example 
given at the beginning of Section 5. We recall again that B* is 
singular and that 


CS (s ) = CS ( s ) 


s+1 1 

_s-2 2_ 
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By Eq. (124) 


Suppose we choose P = 


/N 

C = PC = 

■ 3 

0 

-1 

o" 

- p 

-2 

1 

2 

0_ 



so that P 


-1 



Consequently, PCS (s) = C S (s) 

c 


Note that P cannot alter the zeros associated with |CS(s) 

A /V 

|C S ( s ) | = |pj |CS(s) I , where |p| is, of course, a scalar. 
P 

thus equal to 



i . e . 

B* is 
“P 



a nonsingular matrix. The altered system can therefore be de- 
coupled via linear state variable feedback, and three of the four 
closed loop poles of the decoupled system can be arbitrarily 
specified (Section 5.1). We will choose these poles at s = -1, 
-2, and -3 in order to demonstrate an alternate method ot determ- 
ining F, the feedback matrix. If 

" 0 1 " 

G = B* -1 = 

-P 

-1 3 


we can write the final closed loop transfer matrix of the altered 
system directly; 




(S)B B* 
-m-p 


-1 


1 

s^+3s+2 

0 


0 

1 

s+ 3 


s+4 

s+4 


Since all quantities are known with the exception of 5p(s), we 
solve the above equation for 6p(s); 
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1 1 1 II 1 1 1 1 ■■!!!■■ Hill 111 III I I ii I II 1 1 iiiiii ii »■ 


6„(s) = B B* 
-F ' -xn-p 


s +3s+2, 0 


C p S(s) 


-p ( s ) = 


which yields 


0 1 s +3s+2, 03 -1 


-1 3 0, 


s+3 | s-2 2 


s +s— 6 2s+6 

-F ( s } = 2 

- -6s-24 s +9s+20 


Since 6_(s) = 6 (s) - B FS(s) = 6 (s) - B FS(s), in this example, 

— * h — — m — —m 


and B =1, FS(s) = 6(s) - 6„(s), or 

■“in — — — r 


FS (s) 


-s+9 -2s-6 


5s + 2 3 , -7s-19. 


" 9-1 -6 -2 “I 

.23 5 -19 - 7J 


The closed loop transfer matrix of the altered system can be 
premultiplied by P ^ to yield the closed loop transfer matrix of 
the given system; 


CS ( s ) 6 p 1 ( s ) § m §* 1 


1___ 1 

(s+1) (s+2 s+3 


s+4 
1 s+4 

s+3 


We have therefore succeeded in obtaining a transfer matrix whose 
nonzero entries correspond to those of P \ 

The two algorithms we have now presented do not, of course, 
represent the totality of structure theorem algorithms which can 
be employed to achieve closed loop noninteraction. The 
reader will undoubtedly note that certain basic modifications can 
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be made. For example, two rows of C can be altered instead of one 
and P can be replaced by P(s) in certain cases. The two algorithms 
can also be combined, and the combination modified in a number of 
ways. Space limitations, however, prevent us from covering these 
points in more detail. 

6. A HELICOPTER STABILITY AUGMENTATION SYSTEM 

In this section of the report, we will apply several of the 
results presented earlier in order to design a stability augmen- 
tation system, for a hovering helicopter, based on desired handl- 
ing qualities. In particular, we consider the linearized equa- 
tions of motion (ninth order) of a hovering helicopter, and pro- 
ceed to design a linear, state variable feedback control system 
which satisfies certain conditions. These conditions are express- 
ed in terms of the transfer matrix of the system; i.e., in terms 
of the desired closed loop transfer functions which comprise the 
(6x4) transfer matrix of the system. The steps taken to achieve 
the final feedback design are enumerated here in order to facili- 
tate the review of the design method employed. 

1. A brief discussion of the vehicle which will be 
considered is given and our overall design objective 
is compared with an alternative design. 

2. A corollary to the structure theorem is presented, 
which enables us to deal directly with the physical 
state of the system throughout the entire design pro- 
cess. In this way, our engineering insight can be 
employed at various stages in the design to motivate 
the various mathematical steps taken. 

3. The linearized equations of motion of the vehicle 
are presented, and the physical state of the system 
is defined. 

4. A qualitative discussion of the properties which 

we desire in the final closed loop system is presented 
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and a quantitative interpretation of these properties 
is given in terms of a desired closed loop transfer 
matrix. 

5. The practical considerations which constrain our 
design are presented — first qualitatively and then 
quantitatively . 

6. A design of an initial feedback control system, 
based on the corollary to the structure theorem, is 
presented. Use is made of certain characteristics 
of the particular vehicle considered. 

7. The deficiencies of this initial design are dis- 
cussed and corrected in order to satisfy the con- 
straints imposed. 

8. A final design is then obtained. This design is 
discussed on its own merits, and then compared with 
the feedback system we were hoping to achieve. 

The helicopter we consider is the Sikorsky SH-3D Sea King, 
whose primary mission is to detect submerged submarines through 
the use of a sonar ball, which is lowered into the water while 
the helicopter is hovering at an altitude of approximately 40 ft. 
The linearized equations of motion of this vehicle, as well as an 
alternate feedback stabilization system, are given in reference 
38. We point out that the stabilization system obtained in refer- 
ence 38 is based on linear optimal control, the objective there 
being to maintain the vehicle about a fixed point in the presence 
of disturbances. Our design objective, however, is based on de- 
sired handling qualities; i.e., the ability to design a feedback 
control system which ensures "acceptable manual control" of the 
helicopter active; the equilibrium point (the hovering position) . 

The primary design tool we employ is the structure theorem. As 
we will show, however, certain practical constraints, such as 
limited feedback gains and the inability to combine mechanical 
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inputs, prevent the direct application of the structure theorem. 
Furthermore, when dealing with practical systems, it is desirable 
to avoid any transformation of coordinates (linear combinations 
of state variables) . It is usually preferable to work with the 
actual state variables which define the system. For this reason, 
we will introduce a corollary to the structure theorem. The 
corollary is most useful when a significant number of state vari- 
ables are derivatives of other state variables, a condition which 
occurs quite frequently in practice. In particular, we have the 
following corollary (CIO) : 

CIO: Consider the system x = Ax + Bu , y = Cx. 

Suppose A, an (nxn) matrix, is in multi-input companion 
form as defined by Eqs . (13) , (14) , and (15) ; i.e. , A 

can be partitioned into m diagonal blocks , each a com- 
panion matrix of dimension o^, where i = 1, 2, . . . m. 

Furthermore, suppose that B, an (nxq) matrix, has zero 
e ntries everywhere but the m d^ rows , where d^ = 

Ecr^, K = 1, 2, . . . m. If S(s) is now defined by Eq. (21); 
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a j 

and 6(s) = diag[s x ] - A m §(s) , where A^ is the (mxn) 
matrix consisting of the m-ordered d K ~th rows of A, and 
as the (mxq) matrix consisting of the m-ordered d K ~th 
rows of B, then 


T (s ) = C (si - A) 1 B = CS (s ) <5 1 (s)B m 
Furthermore, if u = Fx + Gw, then 


T_, r (s) = C ( si - A - BF) 1 BG = CS ( s ) 6 1 ( s ) B G 

where 


6 f (s) = 6 (s) - B^FS ( s ) 

Proof : The proof is a direct consequence of theorem Tl and 

corollary C2. 

We now present the linearized equations of motion of the 
SH-3D helicopter. The state variables defining the system have 
been reordered and numbered differently than in reference 38 in 
order to allow direct application of the above corollary. In 


particular , 

we 

consider the 

system x 

— 

Ax + Bu 

' Y 

= Cx , 

where 


-.016 

- 

.05 .0025 

0 

0 

-.0001 

0 

0 

-.0047 


0 

0 

1 

0 

0 

0 

0 

0 

0 


1.97 

0 

-.542 

1 

0 

. 548 

0 

0 

.736 


0 

0 

.00018 

-.3242 

0 

0 

0 

0 

0 

A = 

0 

0 

0 

0 

0 

1 

0 

0 

0 


2.61 

0 

-1.94 

-.163 

0 

-1.96 

0 

.01 

-7.25 


0 

0 

0 

0 

0 

0 

0 

1 

0 


.016 

0 

-.0083 

-.193 

0 

-.0043 

0 

-.303 

5.59 


.0047 

0 

-.0024 

-.0007 . 

05 

-.0025 

0 

.0009 

-.033 
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and 


B 


.05 .005 0 

0 0 0 

- 6.15 .69 0 

0 -.424 0 

0 0 0 

0 - 2.13 21.81 
0 0 0 

0 5.12 .174 

0 .01 .05 


0 

0 

0 

0 

0 


. 3475 

0 

- 7.48 

.022 



0 0 
1 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
1 0 
0 1 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 1 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 1 


The normalized state x and control u are defined as follows: 
x-^ - longitudinal velocity 
x^ - pitch 
x^ - pitch rate 
x 4 - vertical velocity 
Xj. - roll 
Xg - roll rate 
x 7 - yaw 
Xg - yaw rate 
Xg - lateral velocity 

- longitudinal cyclic pitch 
U£ - main rotor collective 
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u g - lateral cyclic pitch 

u^ - tail rotor collective 

Given this open loop system, our objective is to design a 
linear feedback control system based on certain desired handling 
qualities and constrained by certain practical considerations. 

We first discuss the desired handling qualities. 

In particular, we would like to achieve a final "decoupled" 
feedback design where longitudinal velocity and pitch (x^, x 3 , 
and x g ) are affected by longitudinal cyclic pitch (u^) alone. 
Simultaneously, it would be desirable if vertical velocity (x^) 
were affected by main rotor collective input (u^) alone, roll and 
lateral velocity (x 5 , x g , and x g ) were affected by lateral cyclic 
pitch input (u 3 ) alone, and yaw (x^ and x g ) were affected by tail 
rotor collective (u^) alone. In addition to these requirements, 
we would also like to select the poles associated with various 
input/output transfer functions as follows. Longitudinal velo- 
city (y^) , lateral velocity (y g ) , and yaw (y^) should be affected 
by u^ , u g , and u^ via pure integration. Furthermore, the transfer 
functions relating pitch (y ^ ) and roll (y^) to longitudinal cyclic 
pitch (u^) and lateral cyclic pitch (u g ) should represent second 
order systems with critical damping (ref. 32) (£ = .707) and an 

undamped natural frequency (ref. 32) of 3 rad/sec (co n = 3). All 
other specified poles should lie in the half-plane Re s<0. These 
conditions represent the handling qualities which we hope to a- 
chieve, using a linear feedback control system. In terms of the 
closed loop transfer matrix, C(sl - A - BF) ^BG, relating the six 
outputs to the four inputs, these conditions would be satisfied 
if 
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C (si - A - BF) 1 BG = 


where the poles associated with t^ts) lie in the half-plane 
Re s<0 and the numerator polynomials n^^ (s) , r^-^Cs), n^( s )f and 
ng^(s) contain no roots at s = 0. 

As stated earlier, we also have some constraints associated 
with the problem. Actually, there are two constraints, both of 
which are given in reference 38 and repeated here for convenience. 
As stated in reference 38, the amplitude of the controller output 
(feedback gains) is limited to ±10% of the total available range 
of the mechanical input (external input w) . Furthermore, physical 
coupling of the mechanical inputs is not permitted (G must equal 
I) . The first of these two constraints will ensure that the pilot 
can recover successfully from a hardover failure in the augmenta- 
tion system. The second, although not explicitly stated as a con- 
straint in reference 38, ensures a consistent open loop mechanical 
system throughout the entire range of operation of the helicopter. 
The requirement for additional hardware design is also eliminated. 
Mathematically stated, these two constraints imply that u must be 
of the form, u = Fx + w, and all elements comprising the feedback 
matrix F must be less than or equal to 2 in absolute values. This 


n H ( s ) 


s ( s^+4 . 24s+9 ) 


n 21 (s) 


s +4.24s+9 
0 


0 


t 32 ^ 


n 43 (s) 


s +4 . 24s+9 
0 


n 63 (s) 


s (s +4 . 24 s+9 ) 


n 


54 
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latter mathematical constraint is the author's own interpretation 
of the ±10% controller output constraint and is based on the final 
feedback design given in reference 38. 

The problem has now been completely formulated, and our task 
is to design a feedback control system which matches the desired 
closed loop transfer matrix as closely as possible. We note here 
that an exact match in impossible, as we will later show. 

We now recall corollary CIO; i.e. T F (s) = CS (s) 6 ^ (s) B , 

where 6 (s) = S(s) - B FS(s). For this example, CS(s) = I and 
F “ITl “ 

.05 .005 0 0 

-6.15 .69 0 0 

0 -.424 0 0 

0 -2.13 21.81 .3475 

0 5.12 .174 -7.48 

0 .01 .05 .022 

Note that the elements comprising the first and last rows of B m 
are considerably smaller than those contained in the remaining 
rows. Also recall that F, the feedback matrix, is constrained 
in magnitude. Consequently, we conclude that the effect of feed- 
back will be considerably more significant on the middle four 

rows of 6 ( s ) than on the first and last rows. Our design will 

“F 

therefore be based on properly altering these middle four rows 
of <5 (s) . The effect on the other two rows should be relatively 

insignificant. 

Previously, we stated that in view of the constraints, the 
desired transfer matrix could not be achieved via the feedback 
u = Fx + w. We will now be more specific regarding this fact. 

In particular, since CS(s) = I, the closed loop transfer matrix, 
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C(sl - A - BF) 1 B = 6_ 1 (s)B rii . If we now recall the definition of 
B* given in Section 5.1, i.e. 


r 


B* = 



C 


-m 




B 

B 



where f^ = min [ { j rC^A^B ^ 0} , n-1] , we note that for this ex- 
ample, B* = B . Furthermore, B* is invariant under linear state 
- -m - 

variable feedback (refs. 34 and 35). Since G is constrained to 
equal I, every nonzero entry of B* or B^ represents a nonzero 
entry in the closed loop transfer matrix, C(sl - A - BF) of 

the system. However, we can achieve a certain number of zeros in 
the closed loop transfer matrix (partial decoupling) by judicious 
choice of the feedback matrix F. In particular, if we select F 
in such a way that the middle four rows of 6 (s) are diagonal, 

“r 

then seven entries in T„(s) will be zero— namely, those entries 

— r 

corresponding to the zero entries in the middle four rows of B . 

— m 

In particular, suppose F were chosen so that 


A + BF 


V V V 


V V V 


V V V 


0 0 1 

0 -9 -4.24 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


0 0 0 0 

0 0 0 0 

-10 0 0 
0 0 1 0 

0 -9 -4.24 0 

0 0 0 0 

0 0 0 0 


0 0 

0 0 

0 0 

0 0 

0 0 

1 0 

-10 0 


v v v 


v v v 


v v v 


where the v's denote unimportant (for the moment) entries cor- 
responding to the F selected. This (A + BF) is possible to 
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achieve since the middle four rows of constitute a nonsingular 
matrix. The corresponding 6 TT ,(s) would then be 

— r 

fvv v v v v 1 


«F (S) 


0 s 2 +4.24s+9 0 0 

0 0 s+1 0 


0 

0 


00 0 s 2 +4 . 24s+9 0 

00 00 s (s+10 ) 


0 

0 

0 

0 


V V 


V V 


V 


V 


If we now compute the closed loop transfer matrix T_(s): S^fsjB , 

-F -F -m 

it is clear that the (2,3), (2,4), (3,1), (3,3), (374) , (4,1), and 

(5,1) entries of T (s) will be zero. There are two things wrong 

— r 

with this choice of F, however. First, if we actually compute 
F, the feedback matrix corresponding to this choice of 6 (s) , we 

"~r 

would discover that certain entries of F have an absolute value 
greater than 2, a violation of one of the two design constraints. 
Secondly, this choice of 6 (s) would produce a closed loop trans- 

— r 

fer matrix with s(s+10) as - the denominator polynomial of the (5,2), 
(5,3), and (5,4) transfer functions? i.e., the transfer functions 
relating yaw to main rotor collective and lateral cyclic pitch 
would contain pure integrators with relatively high gain. This 
would be an undesirable condition, since the helicopter would yaw 
at constant rate if either of these two controls were displaced 
from their equilibrium positions. At this point in the design, 
the structure theorem proves most useful. We first reduce the 
entries of F by altering the initial requirement that the transfer 
functions relating and y^ to u^ and u^ , respectively, contain 
the denominator polynomial s^ + 4.24s +9. We still maintain 
critical damping (£ = .707), but reduce the undamped natural fre- 
quency to 2.5 instead of 3 rad/sec. This still represents an 
acceptable frequency for the human operator to control. We now 
concentrate on that portion of the 6 (s) block which produces the 
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undesirable integrator transfer functions relating yaw to main 
rotor collective and lateral cyclic pitch. In particular, suppose 
that 


«F (S) 


V V 


V V 


V 


V 


0 

0 

0 

0 

V 


s 2 +3.53s+6.25 00 0 

0 
0 
0 

V V V V 



0 

0 

0 

0 

v 


where a, b, and c are, as yet, unspecified constants. The in- 
verse of the dotted (3x3) submatrix of 6 p (s) is now determined; 


1 

s+1 


0 


0 


\ s . . 1-1 

l$ F (s )| 

i_~ i 

(3,3) -(5,5) 


-c 


s +3.53S+6.25 
-as-b 


s(s+l) (s+10 ) 


s (s+10 ) (s +3.53S+6.25) 


s (s+10 ) 


Postmultiply this expression by the appropriate portion of B m , 
namely 


I o I — 

Ln! 

( 3 ,2 ) - ( 5 , 4 ) 


-.424 0 0 

-2.13 21.81 .3475 

5.12 .174 -7.48 


yields the appropriate block (submatrix) of the closed loop 
transfer matrix, T p (s); 
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, r Vs)| = 


(3 , 2 ) - (5 , 4 ) 


-.424 

s+1 


-2.13 


21.81 


.3475 


s 2 +3 . 53s+6 . 25 s 2 +3 . 53s+6 . 25 s 2 +3 . 53s+6 . 25 


-F 

-5,2 


(s) 


where 


-F 

-5,3 


(s) 


Ip 

-5,4 


(s) 


-F 

-5,2 


(s) = 


(.424c+5.12s+5.12) (s +3.53s+6.25) + 2 . 13 (s+1 ) (as+b) 


s ( s+1 ) ( s+10 ) (s +3.535S+6.25) 


The constant term in the numerator of Tp „ (s) is equal to 2.65c + 

— ~> r 

2.13b +32. If this quantity is set equal to zero we will achieve 
a pole-zero cancellation at s = 0 and, hence, eliminate the pure 
integrator in this transfer function. Similarly, 


T 

-F 

-5,3 


(s) = 


-21.81 (as+b) + .174 (s +3.53S+6.25) 
s (s+10) (s 2 +3 . 53s+6 . 25 ) 


and the integrator can be cancelled by selecting -21.81b + 

1.0875 = 0. Therefore, to avoid pure integration in the transfer 

matrix elements relating yaw to both main rotor collective and 

lateral cyclic pitch, we set b = .0499 and c =-12.115. Further- 

-7 48 

more, if a is set equal to zero, T ^ = sTsTTcT) ' a desirable 

transfer function (includes pure integration) relating yaw to 

tail rotor collective. In summary, thus far we have concentrated 

on the middle four rows of the closed loop transfer matrix 

T,_,(s) = S_^(s)B . We have shown that if 
-F -F -m 
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A + BF = 


v 

0 

0 

0 

0 

0 

0 

0 

v 


v 

0 

- 6.25 

0 

0 

0 

0 

0 

v 


V 

1 

- 3.535 

0 

0 

0 

0 

0 

v 


v 

0 

0 

-1 

0 

0 

0 

12.115 

v 


V 

0 

0 

0 

0 

- 6.25 

0 

-.0499 

v 


v 

0 

0 

0 

1 


0 

0 

v 


v 

0 

0 

0 

0 


- 3.535 0 


0 

0 

v 


V V 

0 0 

0 0 

0 0 

0 0 

0 0 

1 0 

-10 0 

v v 


then the resulting closed loop transfer matrix will possess cer- 
tain desirable features. Before discussing this point any further, 
we will calculate F, the feedback matrix in order to ensure that 
all elements have an absolute value less than or equal to 2 and 
also to evaluate the effect of F on the first and last rows of 
T ( s ) . To obtain F, we note that the third, fourth, sixth, and 
eighth rows of both A and (A + BF) are known. Denoting these 
rows of A, (A + BF) , and B as A , (A + BF) r , and § r , respectively, 
we note that (A + BF) = A + B F; i.e. (A + BF) - A = B F, 

— — X~ “ 2T JL -i 1. x 

and since B r is nonsingualr, F = [ (A + BF) r - A^J . Now, 


(A + BF) r - A r = 


- 1.97 - 6.25 

0 0 
- 2.61 0 
-.016 0 


- 2.993 -1 

-.00018 -.6758 
1.94 .163 

.0083 12.308 


0 -.548 

0 0 

- 6.25 - 1.575 

-.0499 .0043 


0 0 -.736 

0 0 0 
0 -.01 7.25 

0 9.607 - 5.59 


and 
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-.1626 -.2646 

0 

0 





B - 1 = 

-r 

0 

- 2.3585 

0 

0 






0 

-.2045 

.0458 .0021 






0 

- 1.62 

.001 

-. 1336 J 



Therefore , 









‘ . 3203 

1.016 

.4867 

.3414 0 

.0891 0 

0 

. 12 


0 

0 

. 0004 

1.594 0 

0 

0 

0 

0 

F = 










-.1195 

0 

.0888 

.1719 -. 

2863 -. 

0721 0 

-.02 

. 3202 


0005 

0 

.0012 

-.5495 . 

0004 -. 

0011 0 

1.28 

.754 . 

and ’ 

we note 

that all elements of F ; 

satisfy 

the constraint 

that 

|f ij 

1 1 2. 

In fact 

, the 

largest is 

1.594 . 

Once F has been de- 

termined, we 

can compute 

A + BF explicitly; 

i . e . 






0 

0 

.0275 

.017 

0 

.0044 

0 

0 

.0013 



0 

0 

1 

0 

0 

0 

0 

0 

0 



0 

•6.25 

-3.535 

0 

0 

0 

0 

0 

0 



0 

0 

0 

-1 

0 

0 

0 

0 

0 

A + BF 

= 

0 

0 

0 

0 

0 

1 

0 

0 

0 



0 

0 

0 

0 

-6.25 

-3.535 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

1 

0 



0 

0 

0 

12.115 

-.0499 

0 

0 

-10 

0 



-.0013 

0 

.002 

.012 

.04 

-.006 

0 

.0281 

0 

and , by 

inspection 

r 
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s 

-.0275s 

-.017 

-.0044s 

0 

-.0013 

0 

s 2 +3 . 535s+6 . 25 

0 

0 

0 

0 

0 

0 

s+1 

0 

0 

0 

0 

0 

0 

s 2 +3.535s+6.25 

0 

0 

0 

0 

-12.115 

.0499 

s (s+10 ) 

0 

.0013 

- . 002s 

-.012 

. 006s- . 04 

-.0281s 

s 


and 


i ; 1 


(s) 


1 

.0275 

.017 

.0044 


.0013 

S 

s 2 +3 . 535s+6 . 25 

s (s+1) 

s 2 +3 . 535s+6 . 25 

0 

2 

s 

0 

1 

n 

n 



s 2 +3 . 5 35s + 6 . 25 

u 

u 

0 

0 

0 

0 

i 

s+1 

0 

0 

0 

0 

0 

0 

1 

s 2 + 3 . 5 35s+6 . 25 

0 

0 

0 

0 

12.115 

- . 0499 

1 

0 

s ( s+1 ) (s+10 ) 

s (s+10 ) (s 2 +3 . 535s+6 . 25 

s (s+10) 

,0013 

.002 

. 012s+ . 46 

-.006S+.04 

.0281 

1 

2 

s 

s 2 + 3 . 535s+6 . 25 

s (s+1) (s+10 ) 

s (s 2 + 3.535s + 6.25 

s (s+10) 

s 


Once 6 ^(s) has been calculated, the overall (6x4) closed loop 

- -1 

transfer matrix T^Cs) = 6 (s)B can be determined, as follows: 

~~_b — r rn. 


. 0 5s + . 007 7s+ . 3175 
s (s 2 + 3. 535S + 6.25) 

-6 . 15 

s 2 + 3 . 5 35s+6 . 2 5 


-.0123 


s +3 . 5 35s+6 . 2 5 


s +3. 535s+6.2S 


-.424 

s+1 


-2.13 


s +3. 535s+6. 25 


5 . 12s +28 . 36s+68. 37 
(s+1) (s+10) (s 2 + 3.535s+6.25) 


.096 


s + 3 . 5 35s + 6 . 2 5 
0 


s +3.535S+6.25 


.01s +. 297s +1. 05s +.8 04s-. 55 . 0 5s + . 04 6s+ 1 . 19 


s ( s+ 1 ) (s + 10) ( s 2 + 3 . 5 35s + 6 . 25 s (s 2 +3 . 535s+6 . 25) 


0 

. 3475 


s +3. 535S+6.25 


■ 17 4 s+ . 6 1 5 

(s+10) (s 2 +3. 535S+6. 25) 


s (s+10) 


.022S+.01 
s ( s+ 10 ) 


LONG. 

VEL. 


VERT. 

VEL. 


LAT. 

VEL. 


LONG. CYCLIC 
PITCH 


MAIN ROTOR COLLECTIVE 


LATERAL CYCLIC PITCH TAIL ROTOR 
COLLECTIVE 


91 


This T„(s) is proposed as the final closed loop transfer 
matrix. If we compare this transfer matrix with the original 
desired closed loop transfer matrix, we can make several observa- 
tions. First, the requirement that pure integrations appear in 
the transfer functions relating longitudinal velocity to longi- 
tudinal cyclic pitch, lateral velocity to lateral cyclic pitch, 
and yaw to tail rotor collective has been met. Furthermore, by 
appropriate selection of the feedback matrix F, we have elimin- 
ated pure integrations in almost all of the other transfer func- 
tions. The only exceptions are the transfer functions relating 
lateral velocity to the main rotor and tail rotor collective in- 
puts. We note, however, that the gains associated with the in- 
tegrator portion of these two transfer functions are quite low, 
relative to the lateral cyclic pitch input. We were unable to 
achieve the undamped natural frequency of 3 rad/sec relating pitch 
and roll to longitudinal cyclic pitch and lateral cyclic pitch, 
respectively, because of the constraint placed on the magnitude 
of the elements comprising the feedback matrix F. A reduction of 
this desired frequency to 2.5 rad/sec, however, was found to be 
obtainable without excessive feedback gains. Continuing our 
analysis, we note that we have achieved as much decoupling as 
possible, consistent with the constraint that G = I. (This point 
was discussed in considerable detail earlier.) We also note that 
longitudinal cyclic pitch affects pitch and longitudinal velocity 
almost exclusively. Vertical velocity is affected by main rotor 
collective only. Lateral cyclic pitch affects lateral velocity 
and roll, as desired, with some unavoidable cross-coupling into 
yaw. The main input affecting yaw is tail rotor collective, as 
desired . 

Although not all design objectives were met, it appears that 
this design is about as close to the desired design as possible, 
consistent with the constraints on F and G. If there were no 
constraints on F or G, we could have achieved almost a perfect 
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match between desired and actual closed loop transfer matrices. 

The only exceptions would be nonzero entries appearing in all 
elements of the first and last rows of the closed loop transfer 
matrix, in general. 

7. CONCLUDING REMARKS 

We have introduced the structure theorem, a new and rather 
basic idea which was shown to have broad implications in the anal- 
ysis and design of linear multivariable control systems. The 
underlying idea behind the structure theorem is the ability to 
combine time and frequency domain information in a compact and 
concise expression through a transformation (restructuring) of 
the dynamical equations of motion which characterize the system. 

A significant number of new results were obtained, and cer- 
tain prior results were given a new interpretation. In particular, 
in Section 3 we presented a new method for obtaining realizations, 
given the transfer matrix of a linear multivariable system. We 
also derived a direct, constructive technique for arbitrarily 
assigning all n eigenvalues of a controllable multivariable system 
using linear state variable feedback. The question of pole place- 
ment via output feedback was then considered. We presented a 
compact expression for the closed loop characteristic equation 
(closed under linear output feedback) of a multivariable system, 
and then demonstrated how this result could be used to obtain 
certain desired closed loop poles. 

In Section 4 , various questions pertaining to linear optimal 
control were considered. Few new ideas were presented here, but 
a number of improvements were made over previous work. In par- 
ticular, we demonstrated how a solution to the output regulator 
problem could be formulated in the frequency domain. This result 
led to two new expressions for characterizing an optimal feedback 
design. We then demonstrated how spectral factorization could be 
used to obtain a solution to the regulator problem, and concluded 
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by discussing a new technique for formulating the root-square 
locus . 

We considered the design of noninteracting control systems 
in Section 5. A number of new results and design algorithms were 
presented. We gave a procedure for determining what poles could 
not be altered while decoupling via linear state variable feed- 
back. We also presented a design algorithm for achieving maximum 
pole assignment under linear decoupling feedback. We then char- 
acterized those systems which can be decoupled via input dynamics 
and presented an algorithm for achieving an asymptotically stable 
decoupled design with arbitrary pole placement. Two design algo- 
rithms for achieving a noninteractive design via linear state 
variable feedback were then given and demonstrated by example. 

In Section 6, we considered the design of a helicopter 
stability augmentation system based on desired handling qualities 
Many of the results obtained in the earlier sections were applied 
Here, it was shown that pencil and paper methods can be employed 
to achieve an acceptable feedback control system. 

Preliminary computer programs have been devised to implement 
much of the theory developed in this report. Some of these pro- 
grams were mentioned briefly. Currently, there is an effort 
underway to combine these preliminary programs in a highly inter- 
active man-computer system which will facilitate the analysis and 
design of complex multivariable systems. 

Although a significant amount of work has been accomplished, 
much remains to be done. Answers to additional questions are re- 
quired in most of the areas already mentioned as well as in re- 
lated research areas such as model matching, state estimation, 
sensitivity reduction, and nondeterministic systems. The struc- 
ture theorem, which we have introduced in this report, appears to 
offer a valuable design and analysis technique for achieving 
further results in these areas and others. 
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Appendix A 
LIST OF SYMBOLS 


Symbol 

Page 


Definition 

X 

4 

The 

n-dimensional state of a linear system 

A 

4 

The 

constant (nxn) state matrix 

B 

4 

The 

constant (nxm) input matrix 

C 

4 

The 

constant (pxn) output matrix 

Z 

4 

The 

p-dimensional output of a linear system 

u 

4 

The 

m-dimensional input of a linear system 

?o (s) 

4 

The 

open loop transfer matrix 

y_(s) 

4 

The 

Laplace transform of the output y_ 

U (s) 

4 

The 

Laplace transform of the input u 

s 

4 

The 

Laplace operator s = o+jw 

w 

4 

An m-dimensional external input 

F 

4 

A constant (mxn) state variable feedback matrix 

G 

4 

A constant (mxm) feedforward matrix 

5f,g ( s » 

5 

The 

closed loop transfer matrix 

0 

5 

An 

(nxn) nonsingular similarity transformation 



matrix 

z 

5 

The 

altered state under Q, Qx 

A 

5 

The 

altered state matrix under Q, QAQ 

A 

B 

5 

The 

altered input matrix under Q, QB 

<UI 

5 

The 

altered output matrix under Q, CQ ^ 

/\ 

F 

6 

The 

altered feedback matrix under Q, FQ ^ 

A 

m ^ 

“F f G 

6 

The 

altered closed loop transfer matrix under Q 
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Symbol 

Page 

Definition 

K 

8 

The (nxnm) controllability matrix (see K, 
pp. 32 and 58 also) 

L 

9 

The (nxn) nonsingular matrix obtained by a 
"lexicographic" ordering of the first (n) 
independent columns of K 

4 

11 

The (nxn) "companion" form closed loop state 
matrix 

T f (s) 

11 

The closed loop transfer matrix when G = I 

S (s) 

12 

An (nxm) matrix of monic single term polynomials 
in s used in defining the closed loop transfer 
matrix via the structure theorem. 

«F (S) 

12 

An (mxm) matrix of polynomials in s used in 
defining the closed loop transfer matrix via the 
structure theorem. 

B 

-m 

12 

An (mxm) upper-triangular matrix obtained from 
B and used in defining the transfer matrix via 
the structure theorem 

V s > 

13 

The characteristic polynomial of the closed loop 
system 

C*(s) 

13 

Shorthand notation for CS(s) 

N p (s) 

13 

A (pxm) matrix of polynomials in s--the 
numerator of the closed loop transfer matrix 

5 f (s) 

13 

An (mxm) matrix of polynomials in s equal to the 
adjoint of 6p(s) 

/\ 

-m 

14 

An (nxm) matrix obtained from A and used in 
obtaining § F (s) 

-o (s } 

14 

An (mxm) matrix of polynomials in s used in 
defining the open loop transfer matrix via the 
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Symbol 

Page 

Definition 



structure theorem; i.e. 6 (s) = <5„(s) when F = 0 

“O “ r •“ “ 

W 

15 

An r-dimensional subspace of 

*N 

15 

The n-dimensional Euclidean space 

w 1 

15 

The orthogonal complement of W 

fii 

15 

One of the basis vectors of W"*" 

B 

-e 

15 

An extension of B which includes the basis 



vectors, 3^, of W 1 

-e 

15 

An (nxn) similarity transformation matrix 
defined in terms of the extended input matrix 
B 

Vu (s > 

16 

-e 

The characteristic polynomial associated with 
the "uncontrollable" part of the closed loop 
state matrix 

— F , c ^ S ^ 

17 

An (mxm) matrix of polynomials in s associated 
with the "controllable" part of the closed loop 
state matrix 

A 

F 

-e 

18 

An (m+qxn) extension of the -Feedback matrix F 

A 

B 

-e , m+q 

18 

/N 

An (m+q) x (m+q) extension of the input matrix B^ 

§e (s) 

18 

An nx(m+q) extension of S(s) 

T(s) 

20 

A given (pxm) transfer matrix— also shorthand 
notation for T q (s) 

{ -c ' — c ' 

C c >21 

A controllable realization of the transfer matrix 

T* (s) 

21 

A modified version of T(s) used in obtaining a 
controllable realization 

A 

24 

An arbitrary set of (n) complex numbers 

A 

A* 

-m 

25 

A -A /N 

A special form of A m +B m F used for arbitrary 


pole placement via state variable feedback 
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Definition 


l l li 


Symbol 

H 

J 

R 

T 

u* 

K 

CO 

A (s) 

A* (s) 

F*(s) 

6 (s) 

£ (s) 
W(s) 

U 

B* 


8 


P 


V 




Page 

27 An (mxp) constant output feedback matrix 

31 A quadratic performance index 

31 An (mxm) positive definite input weighting 
matrix (see p. 71 also) 

T 

31 The "square root" of R; i.e. R = T T 

31 The optimal control which minimizes J 

32 The positive definite solution to the matrix 
Riccati equation (see K, pp. 8 and 58 also) 

33 A frequency in radians/second 

33 The characteristic polynomial of the open loop 

system 

33 The characteristic polynomial of the closed loop 

optimal system 

35 The optimal state variable feedback matrix 

35 Shorthand notation for A q (s) 

41 A symmetric (mxm) matrix of polynomials in s 

41 A "spectral factorization" of £(s); i.e. 5>(s) = 

W T (-s) W(s) 

T 

41 An orthogonal matrix U U = I 

52 An (mxm) matrix used in decoupling multivariable 
systems 

53 The degree of the polynomial p(s) 

54 The maximum number of specifiable poles under 
decoupling state variable feedback 

54 An (mxm) matrix of polynomials in s obtained by 

factoring C*(s) 
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Symbol 

Page 

Definition 

P(s) 

56 

An (mxm) diagonal matrix of polynomials in s 
obtained by factoring C*(s); i.e. C*(s) = 
P(s)C* I (s) 

2ii ( s) 

56 

The open loop transfer matrix of the system 

• 

x = Ax+Bu , = Cjj£ 

-II 

56 

A modified output matrix C^Q 1 S(s) = C^Cs) 

H 

H 

>4 

56 

The output Cj^x 


56 

A diagonal matrix used in obtaining a decoupling 
state variable feedback matrix for arbitrary 
pole placement 

A* 

56 

An (mxn) matrix used in obtaining a decoupling 
state variable feedback matrix 

G(s) 

57 

The (mxm) transfer matrix of a dynamical system 
used for decoupling via "input dynamics" 

2 

58 

The k-dimensional state of the system whose 
transfer matrix is G(s) 

u 

-e 

58 

The output of the system whose transfer matrix 
is G(s) 

V 

58 

The input to the system whose transfer matrix 
is G(s) 


59 

The state of the composite system consisting of 
G ( s) in series with the given system 

{K,L,M,N} 

59 

The quadruple of constant matrices defining the 
k-th order system whose transfer matrix is 
G(s); i.e. G(s) = M(sI-K) -1 L+N (see K, pp. 8 
and 32 also) 

*c 

59 

The output of the composite system — equivalent 


to y . 
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Symbol 

gggg 

Definition 

Ic (s) 

60 

The transfer matrix of the composite. system 

{ ?C'S C > 

61 

A constant composite feedback pair used to 
decouple and specify poles 

T (s) 

-II 

63 

The transfer matrix associated with a factored 
version of the composite system 

p s (s) 

63 

The "stable" polynomial factor of J C* ( s ) j 

p u (s } 

63 

The "unstable" polynomial factor P u (s) of 
|C* I (s) | ; 1 QxJ ( s) ! =Cp g (s)p u (s) 

B* 

-c 

69 

The matrix B* (see p. 52) associated with the 
composite system 

in > 

tQ 

71 

/\ 

An alteration of the matrix C used in achieving 
noninteraction 

R 

71 

A A A 

An (mxn) matrix added to C to yield C ; i.e. C 

- _ q -q 

C+R (see p. 31 also) 

B* 

-q 

71 

The matrix B* (see p. 52) defined in terms of 

A A 

C rather than C 
-q 

A 

e P 

74 

/N 

An alteration of the matrix C used in achieving 
noninteraction 

P 

74 

/A, 

An (mxm) matrix which premultiplies C to produce 

/\ ^ -A 

C ; i.e. C = PC 

ir ir 

B* 

~P 

74 

The matrix B* (see p. 52) defined in terms of 

/A A 

C rather than C 

“P 

K 

82 

The damping ratio associated with a second order 
linear system 

w n 

82 

The undamped natural frequency associated with a 
second order linear system 
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NOTE: 


Appendix B 
SPECIAL NOTATION 

Underlined letters refer to vectors or matrices 


I 

;P * 

The 

(pxp) identity matrix 

X , 

The 

derivative of x with respect to time 

a. or A. 
-l -l 

The 

i-th row of the matrix A 

a . . 
ID 

The 

ij-th element of the matrix A 

a. . or A. , 
"ID "ID 

The 

ij-th block of a portioned matrix A 

A' or A T 

The 

transpose of the matrix A 

A- 1 

The 

inverse of the matrix A 

a~ t 

The 

-1 . T 

transpose of A or the inverse of A 

adj (A) 

The 

adjoint of the matrix A: i.e. 


A" 1 

= adj (A) v | a| 

det(A) or |a| 

The 

determinant of the matrix A 

diag [a i ] 

A diagonal matrix whose entries are a^ 

[ • ] *A ( s ) 

The 

T 

product A (-s)A(s) 

Re s 

The 

real part of s; i.e. o if s = a+jw 

Im X 

The 

imaginary part of the complex number 


A > 0 

I I 
I A I 
I I _l 

(i , j ) - (k, 1) 


A nonnegative definite matrix A 

The (k-i+1) x (1- j+1) submatrix of A consisting of 

those elements common to the i through k-th 
rows and j through 1-th columns of A 
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